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Abstract: In this article, we solve the loop equations of the /3-random matrix model, 
in a way similar to what was found for the case of hermitian matrices (3 = 1. For 
P = 1, the solution was expressed in terms of algebraic geometry properties of an 
algebraic spectral curve of equation = U{x). For arbitrary jS, the spectral curve is 

no longer algebraic, it is a Schrodinger equation {{hd)"^ — U{x)).ip{x) = where 
h oc (a/5 — In this article, we find a solution of loop equations, which takes 

the same form as the topological recursion found for (3 = 1. This allows to define 
natural generalizations of all algebraic geometry properties, like the notions of genus, 

cycles, forms of 1st, 2nd and 3rd kind, Riemann bilinear identities, and spectral 
invariants Fg, for a quantum spectral curve, i.e. a D-module of the form y"^ — U{x), 
where [y, x] = h. Also, our method allows to enumerate non-oriented discrete surfaces. 
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1 Introduction 

Spectral invariants and algebraic geometry 

In [T2l [T6] , was presented the definition of spectral invariants Fg for any algebraic plane 
curve, i.e. given by a polynomial equation 

£{x,y) = ^Sijx'y^ = 0. 

id 

Those invariants Fg{S) are defined in terms of algebraic geometry quantities defined 
on the Riemann surface of equation S{x,y) = 0. Their definition involves residues at 
branchpoints of some meromorphic forms. Their definition provides a natural basis of 
meromorphic forms of 1st, 2nd and 3rd kind, and a natural framework for all algebraic 
geometry notions. 
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Moreover, the invariants Fg of [T6j have many nice properties, for instance their 
deformations under changes of the complex structure of £ is given by some "special 
geometry" relations, and provide a natural form-cycle duality. Also, they are invariants 
under changes of £ which conserve the symplectic form rfx A in C x C, they have 
nice modular properties, and finally, they define the tau-function of some dispersionfuU 
integrable system associated to £. 

Also, those invariants Fg have deep relationships with enumerative geometry, for 
instance they have been related to the Kodaira-Spencer theory [10], to combinatorics 
of discrete surfaces (maps), to intersection theory [l6l[Tl], and they are conjectured to 
be equal to the Gromov-Witten invariants of some toric Calabi-Yau target 3-folds [3] . 

Algebraic geometry on "quantum" curves 

Here, our goal is to define those notions for a "quantum curve", where £{x,y) is a 
non- commutative polynomial of x and y: 



The notion of quantum curve has arised in many ways in the litterature [9], and is also 
called D-modules, i.e. a space of functions quotiented by Ker y), where y = hd/dx. 
In other words, one has to study functions ip{x) such that: 



In our attempt to define the spectral invariants analogous to those of [TB] for such D- 
modules, we are naturally led to define all analogous properties of algebraic geometry. 
For instance we define the notions of branch points, sheets, genus, cycles, forms, 
Bergman kernel, and so on... 

Because of non-commutativity, some notions like branch-points, cuts and sheets, 
become "blurred" or "non-localized", i.e. the branchpoint is no longer a point, but a 
"region" of the complex plane, and cuts are asymptotic accumulation lines of points. 

But, otherwise, it is surprising to find that almost all relationships of classical 
algebraic geometry, remained unchanged when h 0, for instance the Riemann bilinear 
identity, the Ranch variational formula, and the topological recursion defining the 
spectral invariants. 

Moreover, we shall find, that in order for our quantities to make sense, we must have 
a "vanishing monodromy" condition, which can be interpreted as a Bethe ansatz, and 
this gives a geometrical interpretation of the Bethe ansatz. 





£{x, hdx).ip{x) = 0. 
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Let us also mention that in a previous article [15], we treated a special case, where 
the Schrodinger potential U (x) was quantized, and we shall see, under the light of this 
new work, that it was the case of a degenerate quantum surface, with no branchpoints. 

Hyperelliptical case 

Here, for simplicity, we shall restrict ourselves to polynomials of degree 2 in y (called 
hyperelliptical in algebraic geometry), of the form: 

£ (x, y) = y'^ -U (x) , [y, x] = h 

i.e. to the Schrodinger equation: 

We leave the higher degree case for a further work. 
Link with /? matrix models 

The spectral invariants Fg were first introduced for the solution of loop equations 
arising in the 1-hermitian random matrix model fL2\ [6] . They were later generalized 
to other hermitian multi-matrix models [71 [T7] . 

There exist other matrix models, which are defined with non hermitian matrices. In 
fact it is well known since Wigner [21] that depending on the symmetry of the problem, 
it is sometimes interesting to have matrices that are not hermitian. (For example, real- 
symmetric, unitary, orthogonal or quaternionic, ...). Therefore, it seems reasonable to 
extend the definition of the spectral invariants for those other models. Those other 
matrix models are often called /3-ensembles, and they are classified by an exponent 
(3. The 3 Wigner ensembles (see [21], and we changed (3 l3/2) correspond to /3 = 1 
(hermitian case), (3 = 1/2 (real symmetric case), (3 = 2 (real self-dual quaternion case), 
but it is easy to define a (3 one-matrix model for any other value of (3 (see section M for 
more details). 

In [5], a first attempt to generalize the solution of [12] to other matrix models was 
proposed, but it was not as nice as the topological recursion of [12]. In [S], it was 
assumed that (3 = 0(1) when the size of the random matrix becomes large, and it 
was found that all spectral invariants were related to a double series expansion of the 
form: 

^Ar2-2.-fc (V^-1/V^)'=F,,, 

g,k 

The coefficients Fg^k were computed in [5]. Here, in this article, we shall work at fixed 
h = {^/P — l/^/P)/N, instead of fixed (3, i.e. we shall define the resummed F^'s as: 

Fg{h) =Y, Fg,U. 

k 
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The Fg^kS of [5] can be recovered by computing the semi-classical small h expansion 
of Fg{h). In this article we shall argue that Fg{K) is the natural generalization of the 
symplectic invariants of [16] for a "quantum spectral curve" S{x,y) with [y,x] = h. 

The tool which we use for studying the /3-matrix model, is the loop equation method. 
Loop equations are related to the invariance of an integral under change of variable. 
They can be obtained by integrating by parts. Loop equations for the /5-matrix model 
have been written many times [11], [13], and here we show how to solve them order by 
order in 1/A^, at fixed h. 

The /3-matrix model and its loop equations are explained in section [9l 



2 Schrodinger equation and Bethe ansatz 

2.1 Schrodinger equation, generalities and notation 

Let: 

fftlj"{x) =U{x)^Ij{x) (2.1) 

be a Schrodinger equation with U (x) a polynomial. Let U (x) be a polynomial of degree 
2d, and define the polynomial "potential" V{x) of degree + 1 by its derivative: 

d 

V'{x) = 2 {VU)+ = J2 tk+i (2.2) 

fe=0 

where ()+ means the polynomial part of the Laurent series at x — > oo. We also define: 

P,.) . m _ ui.) - ,2.3) 

so that P is a polynomial of degree d — 1. 
Eventually, we define: 

to = lim ^ (2.4) 

x^oo V [X] 

Remark 2.1 Just in order to give names to those parameters, let us say that in the lan- 
guage of integrable systems, the coefficients to)*i)i2> ■ ■ ■ itd+i are called the "Casimirs", and 
the remaining coefficients of P are the "conserved charges" . They will play a special role later 
on in this article. In matrix model language (see section [9]), ti, . . . are called the times 
associated to the potential V{x), to is often called the temperature, and the remaining coef- 
ficients of P are called "filling fractions". In the language of algebraic geometry, parameters 
tk with k > 1 are coupled to 2nd kind meromorphic differential forms, is coupled to 3rd 
kind, and the remaining coefficients of P are coupled to 1st kind holomorphic differentials, 
see section [6] about form-cycle duality. 
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2.1.1 Stokes Sectors 



From the study of the Schrodinger equation we know that the function iplx) is subject 
to the Stokes phenomenon, i.e. although ip{x) is an entire function, its asymptotics 
look discontinuous near oo. We therefore need to introduce properly the Stokes sectors 
by defining the following quantities: Let 

Oo = Arg(trf+i) 

be the argument of the leading coefficient of V{x). 
We define the Stokes lines going to oo as: 

L.= {./A.gM = --A_ + ,|±i} (2.5) 

Those are the lines where IieV{x) vanishes asymptotically. 
We define the sectors: 

i.e. Sk is the sector between Lk-i and Lk- 

Notice that in even sectors we have asymptotically IieV{x) > and in odd sectors 
we have ReV{x) < 0. 







\ Sector 1 / 






\ / Sector 






Sector 2 \ / 






/ \ Sector 5 






Sector 3 / \ 






/ Sector 4 \ 






Line of accumulation of zeros / \ 


Fi^ 


rure 1: Example 


of sectors for a potential of degree deg\^ = 3, i.e. d = 2. If 


dej 


z,V = d + 1 there 


are 2d + 2 sectors. 
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2.1.2 Stokes phenomenon 

Any solution of a linear equation, is analytical where the coefficients of the equation 
are analytical, and it may possibly have essential singularities where the coefficients 
are singular. Here, U{x) is an entire function with a singularity (a pole), only at oo, 
thus, any solution ip is an entire function with a possible essential singularity at oo. 
The asymptotics of ip near oo are subject to the Stokes phenomenon. This means that, 
although ip is analytical in the whole complex plane, its asymptotics at infinity may 
change from sectors to sectors. 

From the study of the Schrodinger equation it is known that in each sector S^, iplx) 
has a large x expansion: 

Sk X 

and the sign ±, may jump discontinuously from one sector to another as well as the 
numbers A^, Bk,Ck, ■ ■ ■ (and in general, all the coefficients of the series in i at infinity). 

2.2 Decreasing solution 

Let us consider a specific solution iIj{x) of the Schrodinger equation which is expo- 
nentially decreasing in some even sector at infinity. For writing convenience, we will 
choose ip{z) = ipo^z) a decreasing solution in sector Sq. Without further indi- 
cation, ^{z) is now understood to be ipo{z) in the rest of the article. Note that 
this choice is quite arbitrary at the moment, and one should wonder if the quantities we 
are about to compute depend on this choice, but we are presently not able to answer 
this question properly, and leave it for further study. 

An important and useful result is the Stokes theorem which claims that if the 
asymptotics of ip{x) is exponentially small in some sector, then the same asymptotics 
holds in the two adjacent sectors (and therefore ip{x) is exponentially large in those 
two sectors). 

In the general case, (i.e. a generic potential U{x)) our solution tp{x) is decreasing 
only in sector 0, and is exponentially large in all other sectors. But if the Schrodinger 
potential U{x) is non-generic (quantized), then there may exist several sectors in which 
■ip{x) is exponentially small. Due to Stokes theorem, if ip is exponentially small in some 
sectors then it must be exponentially large in the adjacent sectors, this implies that 
there are at most d + 1 sectors in which ip is exponentially small. 

The case studied in [15] was the most degenerate case, such that ip is exponentially 
small in d + 1 sectors. 
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2.2.1 Zeroes of ^fj 

The main difference with our previous article [15] is that we will not restrict ourselves 
to the case where iplx) is a quasi-polynomial which can only be obtained with very 
non-generic potential U (x). Here iplx) is an entire function with an essential singularity 
at oo, and with isolated zeroes labelled sf. 



In particular, the number of zeroes of i/j may be finite or infinite. 

If ip{x) has an infinite number of zeroes, it is known that the zeroes may only 
accumulate near oo, and only along the Stokes half-lines Lj's bordering the sectors 
(see fig|2]). In fact, there is an accumulation of zeroes along the half-line Lj if and only 
if ip is exponentially large on both sides of the half-line. 

For example in the case where ip{x) = ipk{x) is a solution that exponentially de- 
creases in sector k then it implies that there is no accumulation of zeroes along the 
half-lines Lk and Lk-i- 

If U{x) is generic, then ip has an infinite number of zeroes, and the zeroes accumulate 
at oo along all half-lines Lj with j l,2d+ 1 (because remember that is implicitely 
assumed to be tpo which decreases in sector 0), i.e. there are generically 2d half-lines 
of zeroes. The situation is illustrated in fig [21 



Figure 2: The zeroes of ip accumulate near oo along the half-lines bordering sectors 
where ipQ is exponentially large on both sides. In particular, there is no accumulation 
of zeroes along Lq and -^2^+1- 




(2.7) 




'o 



Sector 5 
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If U{x) is non-generic (quantized), then there are additional sectors in which ijj is 
exponentially small, and thus there can be no zeroes accumulating along the two half- 
lines bordering these sectors. Remember that from Stokes theorem, each time we have 
a new sector in which ip is decreasing, we have two half-lines less of zeroes. Therefore, 
the number of half-lines of zeroes is always even, and we call it: 

Definition 2.1 The genus g of the Schrddinger equation is defined by: 

20 + 2 = # half-lines of zeroes (2.8) 

And if ip has a finite number of zeroes (i.e. there is no half-line of zeroes), we define 
g = — 1. We have 

-l<0<rf-l (2.9) 

Note also that the definition of g a priori depends on the choice of the solution %l) = 
ipo since two different solutions of the same Schrodinger equation may have different 
numbers of semi-lines of zero accumulation. 

An exception, is in the special cases g = — 1 where it is easy to see that every choice 
of = ip2k would give the same value of g. 

Indeed, consider Q2k and Q2k' be the genus defined from the solutions ip2k exponen- 
tially small in sector and ip2k' exponentially small in sector 5*2^': 

if 02/c = "1, this means that ip2k is exponentially small in all even sectors, in 
particular in sector S2k', and therefore 'ip2k oc ip2k', and therefore Q2k' = — 1- 

2.2.2 Caseg = -1 

The case g = — 1 was studied in jlSj. This is the case where ip has only a finite number 
of zeroes, it is a quasipolynomial: 

?/'(a;) e^^*^^^ = polynomial. (2.10) 

Notice that in order to diminish g by 1, we need to quantize one parameter of U, 
and therefore to reach g, we need to quantize d — 1 — g parameters. In particular, to 
reach g = —1, we need to quantize d parameters, i.e. P is completely fixed in terms of 
V, and in particular, to is quantized. 

In the applications to random matrices, to is usually a free parameter (called the 
temperature) and is never considered quantized, and therefore the case g = — 1 is never 
obtained in random matrices. 

Another way to say that, is that the case g = — 1 has no h —* classical limit, and 
therefore in classical geometry we always have g > 0. 
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2.3 Resolvent 



The first ingredient of our strategy is to define a resolvent similar to the one in matrix 
models. 

Definition 2.2 We define the resolvent for a generic solution i/j by: 

Jj'ix) V'ix) 
o;(x) = n^ + ^ (2.11) 

It is clear that this function is analytical except at the zeros of ip{x) where it has 
simple poles with residue h: 

h 

u{x) + reg. (2.12) 

It also has a possible essential singularity at infinity with the same location of discon- 
tinuities as ip{x). Eventually, note again that the definition of uj{x) depends on the 
choice of ip{x). 

2.4 Sheets 

In sector Sk we have the asymptotic: 

iIj(x) ^ e^^^""' X s (^^ + + . . .) (2.13) 

Sk X 

where 77^ = ±1. That translates for the resolvent to: 

u{x) ^-±^iV'ix)-h^)-^ + 0{l/x'), (2.14) 

x^ook 2 XX 

Therefore it depends if the solution ip is exponentially big or small in sector k (and of 
course on the parity of k). For a generic g = d — 1 solution which is exponentially big 
in every sector except Sq (and thus has an alternating sign in the exponential) then 
r]k = {-!)'' (except r]o = -1). 

Definition 2.3 We call "physical sheet", the union of sectors where rjk = —1, in 
those sectors we have: 

uj(x) ^- + 0(l/x^) (2.15) 

X 

Notice that the sectors Sq, Si and S2d+i cii^g always in the physical sheet. 

And we call "second sheet", the union of sectors where rjk = +1, in those sectors 
we have: 

uj{x) ^V'{x) + 0{l/x) (2.16) 
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This definition comes from tlie analogy witli the resolvent in matrix model (see section 
[9] for details). 

For a generic potential U{x), all odd sectors are in the physical sheet, and all even 
sectors except 5*0 are in the second sheet. 

Notice that if g = — 1, there is only the physical sheet, i.e. there is no second sheet. 

2.5 The Bethe ansatz 

In the polynomial case studied before [15], a key ingredient for establishing results was 
the Bethe ansatz. This ansatz basically deals with the behaviour of u;(x) around zeroes 
of ip. The zeroes of ip are called "Bethe roots". 

The Bethe ansatz can be formulated in many ways. One way to formulate it, is to 
say that l/ip"^ has no residue at the Sj's: 

in this way, it will play a key role in defining contour integrals, because all integrals of 
the type J dx/ip'^{x) are insensitive to the exact location integration path with respect 
to the Sj's, i.e. such integrals will depend only on the homotopy classes of paths. 

Equation (12.171) can also be formulated, in a form very similar to the Bethe ansatz 
in the Gaudin model [201 [I] cis follows: 

Theorem 2.1 The roots Si of ip satisfy the Bethe ansatz: 

Vi, V'(si) = 2 \im (u(x) —] ■ (2.18) 

x^si \ X-SiJ 

It is a regularized version of the Bethe equation for Gaudin model: 

Vz, V'{s,) " = " 2hy^^— 

when the number of zeros is infinite and the sum is ill-defined, 
proof: 

This theorem is a classical result and is easy, it just consists in rewriting the 
Schrodinger equation as a Ricatti equation. We proceed the same way as in [15] and 
compute: 

V\x)uj{x) - u?{x) - fiuj\x) 



'\l}{x) 2 \ 4 ■\\){x) ip'^{x) 
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V'jxy ^2 ^"{x) ^ V"{x) 
4 ij{x) 2 

= P{x) (2.19) 

which is a polynomial in x, of degree (i — 1. 

From its definition, it is clear that uj"^ + hw' has no double pole at the Sj's, but it 
could have simple poles. Consider now a zero Si of ip, and define: 

CUi{x) = uj{x) 

X Si 

Then, Ui{x) is regular at x = Sj, and we may compute uJi{si). Compute: 



Res (x) + hw'ix) = Res cj. (x) + 2/1—^^ h + ^K^) / \2 

Ui{x) 



= Res 2h- 

= 2htUi{s,) (2.20) 

On the other hand we have, from eq. (I2.19P we have: 

Res u;^(x) + huj\x) = Res V'{x)uj(x) — P{x) 
= ResV'{x)uj{x) 

= hv\s,) (2.21) 

Therefore we find : 

Vi , V'{si) = 2C0i{si). 

This equation is the Bethe equation for the roots Sj's. Note that the potential V'{x) is 
completely determined by the data of the potential U{x) and does not depend on ip. 
In particular, in the case where there are only a finite number of Sj's we recognize the 
Bethe equation for Gaudin model |15j : 



Vz, V'(si) = 2hy" ^ 
which were completely defining the Sj's. □ 



3 Towards a "Quantum Riemann Surface" 

From the definition of our non- commutative spectral curve (i.e the Schrodinger equa- 
tion), it is tempting to generalize the classical notions kwown in algebraic geometry 
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and Riemann surfaces to our "quantum" case ("quantum" is not to be understood as 
"quantized" but as "non-commutative" = h). For a Riemann surface, the central 
notions are those of cuts, sheets, genus, cycles and meromorphic differentials forms of 
1st, 2nd and 3rd kind. In our context, the picture needs a proper adaptation in order 
to recover the terminology of Riemann surfaces and algebraic geometry. 

In this section we will define the notions of genus, v4-cycles, S-cycles and the first 
kind differentials dual to them. Here, let us assume that g > 0. 

3.1 Cuts 

First, we like to think of the 2 sheets, as the sectors which correspond to the 2 possible 
behaviors of the resolvent at oo: uj{x) ~ to/x (physical sheet) or uj{x) ~ V'{x) (second 
sheet). 

Then, we consider the cuts as sets of roots Sj's. In some sense, each pair of half 
lines of accumulation of zeroes can be thought of as a cut. 

Definition 3.1 We define cuts as pairs of half-lines of zeroes. 

There is some arbitrariness in grouping the half-lines of zeroes by pairs. 

There is g + 1 cuts, like in classical algebraic geometry, and notice that the case = — 1 
which has no classical counterpart, has no cuts. 

Notice that, contrarily to classical geometry, where the endpoints of the cuts are 
zeroes of U{x), here the endpoints are somehow blurred, we may move a finite number 
of Si's from one cut to another. 

3.2 Cycles 

In standard algebraic geometry, the non-contractible ^-cycles are often thought of as 
surrounding cuts in the physical sheet, and their dual fi-cycles are going through the 
cuts, from one sheet to the other, see fig [31 

3.2.1 A-Cycles 

Consider the complex plane from which we remove the second sheet (sectors where 
uj{x) ~ V'{x)). It is clear that it contains + 1 sectors near oo, and there are g 
homologically linearly independent contours which link them. 

Definition 3.2 We define A-cycles Ai,...,Aq as g linearly independent non- 
contractible contours going from oo to oo in the physical sheet. 
A choice of A-cycles is not unique. 
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— ^ A cycle 












Non independant A cycle 
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' ' y''^ B cycle 
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Representation in two siieets of a Riemann surface of genus 1 


Figure 3: Representation in 


two sheets of a Riemann surface of genus 1. 



Remark that this notion really makes sense only for g > 1. 

Notice that each time i\}{x') ~ g-^(^)/2'i. even sector, it means it is exponentially 
small and thus it also behaves like e~^*^^)/^'^ in the neighboring odd sectors. That means 
that we can always choose v4.-cycles going from odd sector to odd sector. 

Since the first sheet and second sheet are separated by half-lines of accumulations 
of zeroes, every ^-cycle surrounds an even number of such half-lines of accumulations 
of zeroes, i.e. surrounds the cuts in the physical sheet. Like in standard algebraic 
geometry, the cuts are identified as pairs of half-lines of zeroes accumulations and the 
A cycles are going enclosing these cuts. 

3.2.2 Examples 

In the generic case g = c/ — 1, we can define d ^-cycles but only d — \ are linearly 
independent. See picture where d = 7: 
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exp{+V/(2h)) 

We clearly see that the dashed contour is not linearly independent with the others 
since the global sum of the contours (dashed included) is contractible in the physical 
sheet. 

For a non-generic case, there are sectors at infinity where ip is exponentially small. 
In these cases, the definition of the contours need some adaptations because these 
sectors correspond to "degenerate" cuts. Here are a few examples of how to deal with 
these cases. Basically, each time there are two sectors where is small we can replace 
one of the standard A cycle, by a ^ cycle (sometimes called also " degenerate" A cycles) 
that connect them. Here are some examples of the contours in more and more peculiar 
situations for d = 7: 
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exp{+V/(2h)) 



From then it is easy to generalize into more complicated frames: 




exp(+V/(2h)) 



15 



It is then easy to generalize the method in more sophiticated situations. 

In the extreme case where ip is exponentially small in all even sectors, there are 
only d independant "degenerate" A cycles and no A cycles, the genus is = — 1. This 
is the polynomial case studied in [13] where there are no A cycles. 

From the definitions, it is easy to see that the genus q defined above corresponds 
to the number of independant A cycles (we exclude the A cycles). It is also obvious 
that the sum of independant A and A cycles always equals d — 1. 

3.2.3 B-Cycles 

As in classical algebraic geometry, it is standard to define the B cycles with an origin 
lying in the non- independant cut. Moreover, although it would be possible to define 
B cycles attached to the A cycles, we prefer limiting ourselves to the definition of B 
cycles attached only to the A cycles. Basically, they start from the non-independant 
cut, goes through their corresponding A cycle and end at infinity in the same sector 
as their corresponding A cycle. As there are two sectors in which their corresponding 
A cycle ends, we double them so that one goes into one sector and the other one in 
the second sector. We also choose the whole so that they intersect only with their 
corresponding ^-cycles: 

AanB/3 = 26a,f, (3.1) 

This definition is easier understandable with the following pictures: 
Generic case: 
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cxp(+V/(2h)) 



And in a degenerate case: 




exp(+V/(2h)) 
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3.3 First kind functions 

After defining the cycles, another important step is to define the equivalent of the first, 
second and third kind differentials. In this section, we propose a definition of the first 
kind differentials. 

Let hk,k = l,...,d — l,hea basis (arbitrary for the moment, but we will choose it 
orthonormal later on), of the complex vector space of polynomials of degree < d — 2. 
To have more convenient notation, we will label the ^-cycles as Aa ,0+l<«<(i— 1 
and the standard A are labelled Aa, 1 < a < 0. 

Consider the following functions: 

Vk{x) = —-r—- / hk{x')ijj^{x') dx' , deghk<d-2. (3.2) 

Notice that, thanks to the Bethe ansatz, Vk{x) has double poles with vanishing 
residues at the s/s (the zeroes of ip), and behaves like 0(l/x^) in sector So and in 
sectors where is exponentially large, (because the polynomial is of degree less than 
d — 2). Therefore, the following integrals are well defined: 



h,a= f Vk{x)dx , a = l,...,g, /c = 1, 1. (3.3) 

J Ac 

For the degenerate contours Aa, we cannot take the integral since it would not 
converge. We define instead: 

Ik,a= / hk{x) il)^ [x) dx , a = Q + 1, . . . ,d — 1, k = 1, . . . ,d — 1. (3.4) 

J Aa 

The matrix „ with A;, a = 1, . . . , c? — 1 is a square matrix, which gives a pairing 
between the set of paths { Aa,Aa } and the space of polynomials of degree at most 
d — 2. Let us choose a basis hk, dual to the ^-cycles, i.e.: 



(3.5) 



Choosing this set of polynomials gives then the following relations: 

W i = 1, . . . , Q, j = 1, . . . , d — I , (f Vj{x) dx = 5ij (3.6) 

(3.7) 



A^ 



yi = Q + I, ■ ■ ■ , d — 1, j = 1, . . . , d — 1 : / hj{x) ip'^{x) dx = 5j^i 

Ja, 



Moreover, from the definitions, we get an asymptotic expression of Vk{x) at infinity: 
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Theorem 3.1 The functions Vk{x) with k < g are such that: 

k=l,...,Q, Vkix) = 0{x-') (3.8) 

in all sectors at infinity. 

And the functions Vk{x) with Q + l<k<d — 1 are such that: 

k = Q + l,...,d-l, Vk{x)=0{x~^) (3.9) 

in all sectors except in the sector where Ak ends, where we have: 

v,{x) = —^ + 0{l/x'). (3.10) 

proof: 

In sector ooq, we clearly have Vk{x) ~ 0{x'^'^^^''^'^) = 0{x^^). And in a sector Si 
where ijj is exponentially small we have: 



Vk{x) 



h^ij^{x) 

and due to our choice of basis eq. (13.51) , we have 



hk{x')ij^{x')dx' + I hk{x')^^{x')dx' 



(3.11) 



in sector Si. □ 

We claim that the function Vk{x) /c = 1, . . . , g are the generalization of holomorphic 
forms (1st kind differentials). 

Remark 3.1 Classical limit. 

The small h BKW expansion -0 ~ e^as I ^ gives: 

±hk{x) 

mx) == (3.13) 

and Vk{x)dx are indeed the holomorphic forms on the algebraic curve = U{x). 

3.4 Riemann matrix of periods 

An interesting quantity in standard algebraic geometry is the Riemann matrix of pe- 
riods which is the integrations of the holomorphic differentials over i3-cycles. Now 
that we have defined properly the cycles, we can define a similar "quantum" Riemann 
period matrix Tj^, z, j = 1, . . . , g by: 

Tjj =^ ® Vj[x) dx. (3-14) 
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Note that this definition makes sense since Vj{x) (j = 1, . . . , g) behaves as 0(l/x^) 
in the sectors where the ^B-cycles go. Also, thanks to the Bethe ansatz, vj has no 
residue at the roots Sj's, therefore those integrals depend only on the homology class 
of S-cycles, and not on a representent. 

Like for the classical Riemann matrix of periods we have the following property: 

Theorem 3.2 The period matrix r is symmetric: Tij = Tj^i. 
proof: 

We anticipate on results which shall be proved later, but which don't depend on 
this theorem. The proof comes directly from theorem 14.101 below, since: 



(h dx (p B(x, z)dz = 2iTx (p dxva{x) = 2i'KTp^oi 
and from the symmetry theorem 14.111 for the Bergman kernel B{x,z) = B{z,x) 



(b dx (p B(x,z)dz = f dz d) dxB{x,z) = 2i'K (p 

JBr JBry Jb„ JBr Jb 



X, z) = 2i7l (h dzVaiz) = 2i7lTn Ft- 



Be 



□ 



3.5 Filling fractions 

In random matrices, the notion of filling fractions, is just the ^-cycle integrals of the 
resolvent. Here, we easily generalize it by the definition: 

Definition 3.3 The filling fractions ei, . . . , are defined as follows: 

« = 1, . . . , g, = ^ £ - ^) + (3-15) 



where the integer Ua is half the number of Stokes half-lines surrounded by the cycle A 

2d+2 



In other words, corresponds to the angular fraction of the complex plane defined 



by the cycle Aa- 

For a = g + l,...,(i— 1 we define 

a = g + 1, . . . , c? — 1, = (3.16) 

And for a = d, we choose a non-independent A-cycle Ad, which surrounds all the Si 's 
which are not surrounded by Ai, . . . , Ag, and define: 
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Note that this definition makes sense because all the cycles Aa go from an infinity 
where u;{x) — ^ ~ O (^). Note also that this definition depends on the exact locus 
of the contour Aa and not only on its homotopy class, since uj{x) has simple poles at 
the Sj's with residue h. If we deform the contour Aa, the filling fractions can change 
by some integer times h. 

In other words, the filling fractions are "blurred" when h ^ 0, they are defined 
modulo an integer times h. In the classical limit ^ — ^ 0, they become deterministic. 

We have: 
Theorem 3.3 

d 

^ea = h (3.18) 

proof: 

When we perform the sum over the contours Aa-, the contour Ad was defined as the 
"complementary" of the others, i.e. so that the sum is contractible. Since the function 
X uj{x) — to/x is integrable at infinity, we find that its global integral is null. With 
the same argument, it is easy to see that Ylt=i = {d+ 1) because we take all Stokes 
lines once and only once. Therefore we get: 

d t 

a=l a=l 

Note that it also tells us that only c? — 1 of the epsilon's are independant. □ 

Remark 3.2 In the case g = — 1, the only filling fraction is ed = to, and it is also the sum 
of residues of uj at the Sj's: 

(^d = to = y2 = h #{si} 

i 

This shows again, that g = — 1 corresponds to a case where io is quantized, namely to is an 
integer times h: 

to/h E N. 

4 Kernels 

One of the key geometric objects in [15] and in [16], is the "recursion kernel" K{x, z). 
It was used in the context of matrix models, to find a solution of loop equations. Here, 
it will also allow us to define the 3rd and 2nd kind differentials. 
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4.1 The recursion kernel K 

First we define: 

and for eacli a = 1, . . . , 0, we clioose a point Pa G Aa and we define: 

where in the last integral, the integration contour between and x", is along Aa- 
This is described in figJH 



A_alpha cycle 




X in A_alpha 

Infinity sector 

Figure 4: Picture of the path of integration used for the definition of the kernel K{x, z) 



For each a = g + l,...,(i— 1, we define: 

C dr' 
Ca{z)= / ^'(x')-F^. (4.3) 

We now need to describe the domain of definition of these functions. 
First, one can see that for a fixed x, these functions are defined for z outside of 
some " cuts" (see figure H]) 

• Choose a path between ooq and x, then K{x, z) is defined for z outside of this 
path. Across the path ]ooo,a;], K{x, z) has a discontinuity: 

• For each a = 1, . . . g, choose a path between ooq and P^-, then Ca{z) is defined 
for z outside of this path, and outside Aa- Across the path jooo^-fo], Ca{z) has a 
discontinuity: 

,C„W = ^— £^ (4.5) 
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and across the path Aa, Ca{z) has a discontinuity: 

SCJ.) = j^^ (4.6) 

• For each a = g + 1, . . . ,d — 1, Ca{z) is defined for z outside of the path Aa- 
Across the path Aa, Ca{z) has a discontinuity: 

SCa{z) = 2i7r tlj\z) (4.7) 

From these remarks, we now define the recursion kernel K{x, z) by: 

Definition 4.1 Definition of the recursion kernel: 

d-l 



K{x, z) = k{x, z)-J2 M^) Ca{z) {U 



a=l 

it is defined for z outside the cuts mentionned above. 

For a fixed z, the analytical properties in x of K{x, z) are the same as those of 
K{x, z) since all Va{x) are analytic. For a fixed z, the primitive of ip'^{x') can be 
defined locally but not globaly on the complex plane. In fact there is a logarithmic cut 
to be arbitrarily chosen on ]ooo, z]. Anywhere out of this cut the function x — > K{x, z) 
is analytic. 

4.1.1 Properties of kernel K 

The definition of the kernel K{x, z) might seem arbitrary at first glance. But in fact, 
the main reason for the introduction of such kernel is that it has many interesting 
properties: 

It is clear from our definitions that: 

Theorem 4.1 For a given z, the kernel K behaves like: 

K{x,z)^ 0{x-^) (4.9) 

when X — > oo m all sectors. 
proof: 

The result is obvious for sector So and for sectors where ip is exponentially big. 
When it is not, the fact that we substract Ca,o: = Q + I, ■ ■ ■ ,d — 1 gives the result. □ 
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Theorem 4.2 We have in all sectors at infinity : 

K{x,z) ~ 0{z^'^). (4.10) 

More precisely we have: 
with 



k=d-l 



Kk{x) = r x'^ ij\x') dx', (4.12) 



ooo 



and 

P d—l „ 

Kk{x) = Kk{x) - Va{x) j) kk{x) dx - ^ Va{x) j> 

proof: 

It is clear that 

" Kk(x) 



iIj\x') x'^ dx' . (4.13) 



fe=0 

where 



ku{x) = -^— I x"'^\x')dx\ (4.15) 



ooo 



hi/j'^{x 

and therefore 

p d—l p 

Kk{x) = kk{x) f kk{x)dx- Vo,{x) (b ij\x')x"'dx' (4.16) 

Now, ii k < d — 2, notice that x'^ is a polynomial of degree < d — 2, and it is thus a 
linear combinations of /ia(a^)'s: 

d-i 

x" = J2bk,phf,{x') (4.17) 

/3=1 

This implies: 

kk{x) = Y,bk,pM^) (4-18) 

13=1 

Taking now the integral over an A cycle and using the normalization choice of hk{x) 
gives: If a < 

/ kk{x') dx' = h,a (4.19) 

and if a > g 

^\x')x"'dx' = bk,a (4.20) 

Ac 
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This implies that Kk{x) = ii k < d — 2, and therefore 

K{x,z) =0{z-'^). (4.21) 

□ 

Theorem 4.3 Let a = 1, . . . , g, and z on the side of Aa which does not contain ooq, 
then: 

i K{x, z)dx = (4.22) 

J Ac 

proof: 

Notice that if z is on that side of Aa, we have Ca(z) = <f. K(x, z) dx, and therefore 
K[x, z) dx = 0. In fact one can see that the addition of the part with the Ca{z) 
was just put there to cancel out the ^-cycle integrals. □ 

4.2 Third kind differential: kernel G{x, z) 

The second important kernel to define is the equivalent of the third kind differential. 
In [15] this kernel was computed from K by derivation, and we use the same definition. 



Definition 4.2 We define the kernel G{x, z) by: 



,2,.,,, K{x,z) ^^ij'iz) 



G{x, z) = -nV{z) = 2h^K{x, z) - hd, K{x, z) (4.23) 

From an easy integration by parts we find: 

^ 1 ^ 2 f dx' , ff{x') ^'(Z) 



x — z i/j'^{x) J^^x' — z \iIj{x') 'ipiz) 

~hY^Va{x)n^)d.^ 

(4.24) 

which shows that near x = z we have G{x, z) ~ ^, i.e. there is a simple pole of 

residue 1 at z = x. Note in particular that -^j^ { t(){x') ~ singularity at 

x' = z and therefore for a fixed z, there is no more any logarithmic cut ]oo, z] as we 
had for K{x, z). 

Note again that a priori, this function of z has the same lines of discontinuity as 
the kernel K{x, z). But notice that the definition of G ensures that all discontinuities 
of K which are proportional to ip^{z) cancel. 
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Theorem 4.4 G{x, z) is an analytical function of x, with a simple pole at x = z with 
residue —1, and double poles at the Sj's (zeros of ip{x)) with vanishing residue, and 
possibly an essential singularity around oo. 

G{x,z) is an analytical function of z, with a simple pole at z = x with residue +1, 
simple poles at z = Sj, and with a discontinuity across Aa-cycles with a = 1, . . . ,Q 
(and thus no discontinuity on Aa): 

SG{x,z) = -2i7rva{x) (4.25) 

proof: 

K{x,z) is discontinuous when z crosses either ]cx3o,a;], ]cxDo,-Pa] or Aa- However, 
the discontinuity of K{x,z) across ]cxDo,a;], ]ooo,-Pa], and Aa is proportional to ip'^{z), 
and this means by derivation that G{x, z) is not discontinuous there. Across Aa with 
« < 0, the discontinuity of K{x, z) is given by eq. (14.71) . and thus, the discontinuity of 
G{x,z) is 5G{x,z) = —2iTxVa{x). 

Since K{x, z) is regular when z = sj, then it is clear that G{x, z) has simple poles 
at z = Sj, with residue —2hK{x, Sj). 

In the variable x, it is clear from the definition and from the Bethe ansatz 12.181 
that K{x, z) has double poles at x = Sj without residue, and this properties follows for 
G{x, z). 

□ 

Theorem 4.5 

G{x,z) = 0{l/x^) (4.26) 

when X ^ oo in all sectors. 
And at large z in sector Sk: 

lim G{x, z) = G{x, oofc) = r]^ td+i Kd-i{x) (4.27) 

2->OOfc 

where r]k = ±1 is such that ip ~ QnkV/2h sector Sk- 
proof: 

The large x behavior follows from theorem 14.11 The large z behavior is given by 
theorem |4.2[ i.e. G{x,z) ~ rikV'{z)K{x, z) ~ rjktd+iKd~i{x)- The sign depends on the 
behavior of the solution in this sector. □. 

Theorem 4.6 Let a = 1, . . . , g, and z on the side of Aa which does not contain ooq, 
then: 

f G{x,z)dx = (4.28) 

J Ac 

proof: 

Immediate from theorem 14.31 □ 
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4.2.1 Semi-classical limit 



We claim that this kernel is the quantum version of the third kind differential. Indeed, 
in classical algebraic geometry a third kind differential is characterized by analyticity 
except a simple pole with non vanishing residue and a proper normalization on A- 
cycles. Here, apart from the discontinuity along the v4-cycles which is expected since 
these contours represent the "quantum cuts", we have analyticity (apart from the Sj's 
which also define the cuts), a simple pole with residue and a good normalization on 
^-cycles. 

In the BKW semiclassical expansion we have iJj ~ e^^^^ and thus 

2 1 

K(x,z) ; (4.29) 

and 

K{x, .) ~ -i- -^== - Mx) C^{z) (4.30) 

and 



G(x, z) ~ 2v^i^(x, z) ~ ^ - 2 v^{x) C^{z) VTT^ 

X z ^[J[X) ^ 



(4.31) 



The form G(x, z^dx has thus a simple pole at x = 2, in the physical sheet with 
residue +1 and in the other sheet with residue —1, and it is normalized on ^-cycles 
£^ G{x, z)dx = 0. This is indeed the usual 3rd kind differential in classical algebraic 
geometry. 

4.3 The Bergman kernel B{x,z) 

In classical algebraic geometry, the Bergman kernel is the fundamental second kind 
differential, it is the derivative of the 3rd kind differential, and it is another major tool 
in classical algebraic geometry. Following the same definition as in [15], we define: 

B{x,z) = -^d,G{x,z). (4.32) 

The kernel B is going to be called the "quantum" Bergman kernel. 

Theorem 4.7 B{x,z) is an analytical function of x, with a double pole at x = z 
with no residue, and double poles at the Sj 's with vanishing residues, and possibly an 
essential singularity around oo. 

B{x,z) is an analytical function of z, with a double pole at z = x with no residue, 
and double poles at the Sj 's with vanishing residues, and possibly an essential singularity 
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around oo. In particular it has no discontinuity along the A cycles, it is 
defined analytically in the whole complex plane except at those double 
poles. 

proof: 

Those properties follow easily from those of G{x, z) of theorem 14.41 In particular, 
it is important to notice that the only discontinuity of G{x, z) is along the ^-cycles, 
and is independent of therefore B{x, z) has no discontinuity there. □ 

4.3.1 Properties of the Bergman kernel 
Theorem 4.8 

B{x,z) =0{l/x^) (4.33) 

when X ^ oo in all sectors. 
And 

B{x,z) = 0{l/z^) (4.34) 

when z ^ oo in all sectors. 
proof: 

Follows from the large x and z behaviors of G(x, z). □ 
Theorem 4.9 B satisfies the loop equations: 

where p!f'\x^ z) is a polynomial in x of degree at most d — 2. And 

(2||| + 9.) [Bix,z)-^^-^yd. =Pf(.,x) (4.36) 

where p2^\z, x) is a polynomial in z of degree at most d — 2. 
proof: 

This theorem is crucial for all what follows, and its proof is rather non-trivial. Since 
it is very long and technical, we present the proof in appendix \M Those equations are 
indeed the loop equations for the 2-point function in the j3 matrix model, see section 

m □ 
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Theorem 4.10 We have for every a = 1, . . . , g; 



and 



B{x, z)dx = , (b B{x, z)dz = {} (4.37) 

J) B{x, z) dz = 2mvaix) (4.38) 
proof: 

The vanishing of ^-cycle integrals in the x variable is by construction and can be 
seen as the consequence of the same result known for G{x, z) on one side of A and 
the fact that B{x, z) has no discontinuity along the ^-cycles. (Therefore, the nullity 
extend on both sides which no longer need to be treated separately). 

For the z variable, notice that if Aa =]oOj,oOj[ goes from ooj to ooj, where both 
oOj and oOj are in the physical sheet, we have: 

I B{x,z)dz= j ' B{x,z)dz = -]-{G{x,oOj)-G{x,oo.-^) (4.39) 

J Aa JoOi ^ 

and from theorem l4.5l (j(a:. oOj) = r]itd+\Kd,-\{x), we get: 

I B{x,z)dz= I ' B{x,z)dz = '^^^ta+iKd-i{x) (4.40) 

J Aa JoOi 2 

and since oOj and oOj are both in the physical sheet we have rji = rjj = —1, and therefore 

B{x,z)dz = 0. (4.41) 

Aa 

And similarly, when performing the integral over Ba, the contribution from infinities 
cancels out since the contour goes in the same sheet. But since Ba intersects its 
corresponding Aa (and only this one) where the primitive —■^G{x,z) is discontinous, 
the result is the jump of G{x, z) along this Aa, that is to say mva{x). Eventually, since 
Ba and Aa intersect twice, we find eq. (14.101) . □ 

One of our key theorems is: 
Theorem 4.11 B{x,z) is symmetric 

B{x,z) = B{z,x) (4.42) 

proof: 

The proof relies essentially on the fact that B{x,z) satisfies the loop equation in 
the two variables. We have: 



ijj{z) %jj{x) 2(x — z] 



2 ' 
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(4.43) 





i,'{z) 




i,{z) \ 


X 


— Z ) 




i,'{z) 






X 


— z J 



This implies: 



xp[z) X — z 

ip[X) X — z 

^ U{x)-U{z) U'{x) + U'{z) 

{x — zY X — z 



(4.44) 



and therefore: 



{x - zf (2^ + d,) Pf (x, z) + 2U{z) + (x - z)f/'(2) 

= (x - ^)2 (2^ + 9,) (2,x) + 2[/(x) + (^ - x)f/'(x) 
ip{x) 

= R{x,z) (4.45) 

Here, the first hne is a polynomial in x, whereas the second line is also a polynomial 
in z. Therefore, R{x, z) is a polynomial in both variables, of degree at most d in each 
variable. Moreover, we must have: 

R{x,x) = 2U{x) (4.46) 

Therefore we must have: 

R[x, z) = ^(- V\x)V\z) - h ^'^"^^ ~ ^'^^'^ - P(x) - P{z) \ + (x - zYR{x, z) 

(4.47) 

where P(x, z) is a polynomial of both variables of degree at most c? — 2 in each variable. 
Putting this back into 14.451 and using the symmetry x ^ z it implies that: 

(2^ + d,) (Pf (x, z) - Pf (x, z)) = P(x, z) - R{z, x) (4.48) 

Then, we can decompose the r.h.s into the basis ha{x)h/3{z) introduced in 13.61 

P(x, z) - R{z, x)= J2 (^--/3 - Rp,a)ha{x)hp{z) (4.49) 
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Integrating the differential equation eq. (14.481) then gives: 

Pf (x, z) - Pf (x, z)=Y, (Kp - Rp,e.)K{x) vp{z) + Ai{x) (4.50) 

where Ai{x) is some integration constant. 

Then using the loop equations 14.91 we find by substraction that: 

and again, integrating this differential equation we find: 

d-l 

P(x, z) - B{z, x)= J2 - Rp,a)va{x) Vfs{z) + A{x) + A{z) (4.52) 



dy ] {B{y, z) - P(z, y)) = p!f\y, z) - P^'\y, z) (4.51) 



where {211)' /ip + d)A = Ai, and A{z) is some other integration constant. 

The large x and large z behavior of B imply that A{x) = A{z) = 0. We thus get: 

B{X, Z) - B{Z, X) = ^(Pa,/3 - Rf3,a) Va{x)vp{z) (4.53) 

Then, using theorem 14.101 



P(x, z)dx = 0= f B{x, z)dz (4.54) 
We find: 

Va,/3, Ra,f3 = R(3,a (4.55) 

that is to say by 14.531 that the Bergman kernel is symmetric. □ 

We claim that all these properties are essential to name this function a "quantum 
Bergman kernel" . Indeed, the symmetry is absolutely necessary and is completely non- 
trivial. The fact that B{x, z) has no discontinuity is also essential since in standard 
algebraic geometry, it is defined everywhere on the Riemann surface. Using all these 
kernels and their properties, we can then generalize easily the recursion of [121 fT5] 
defining the correlation functions. 

4.4 Meromorphic forms and properties 
4.4.1 Definition of meromorphic forms 

Definition 4.3 A meromorphic form TZ{x) is defined as: 

n{x) = —r— / r{x')tp'^{x')dx' (4.56) 



hip'^{x) 



ooo 
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where r(x) is a rational function of x, which behaves at most like 0{x'^ ^) at large x, 
and whose poles rj are such that: 

Res ilj^{x)r{x) =0 (4.57) 

and for all degenerate Aa cycles 

I i;^{x')r{x')dx' = 0. (4.58) 

J Ac 

It is easy to see, that with this definition, the holomorphic forms Va{x), the kernels 
G{x, z) and B{ meromorphic forms of x. 

4.4.2 Analiticity properties 

A meromorphic forms TZ{x), has poles at x = the poles of r(x), with degree 1 less 
than that of r, it behaves like 0(x~^) in all sectors of the physical sheet. From the 
Bethe ansatz, it has double poles at the Sj's, with vanishing residues. 

In particular, it has an accumulation of poles along the half-lines Lj of accumulations 
of zeroes of ip. 

Also, notice that the following integrals are well defined, and independent of homo- 
topic deformations of Aa (in particular independent of where are the Sj's): 

f n{x)dx. (4.59) 

JAc 

4.4.3 The integration contours around branch-points 

Let us choose some contour Ci, such that each Ci surrounds (in the trigonometric 
direction) a half-line Lj of accumulation of zeroes. In other words it surrounds a 
"branch point". Let us also assume that surrounds all roots of ip, i.e. each root 

of ip is enclosed in one Cj. We also assume that contours Ct and Aa do not intersect 
(they have vanishing intersection numbers): 

Vi = 1, ...,20 + 2, ,Va = 1, Cir\Aa = (4.60) 
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exp(+V/(2h)) 



II 

exp(-V/(2h)) 



12 

exp(+V/(2h)) 



4.4.4 Riemann bilinear identity 

For the Riemann bilinear identity, we need the following useful lemma, which we shall 
use very often in this article: 

Lemma 4.1 For every analytical function f{x) which behaves at infinity at most like 



contour Ci ( and thus must he regular at the root Sj 's ) we have, for xq outside of all 
A-cycles (i.e. on the same side as ooq) : 



proof: 

Clearly, the contours Ci enclose no singularity of K{xq, x)f{x) and can be contracted 
to 0. □ 

Then we can write the bilinear Riemann identity: 

Theorem 4.12 Riemann bilinear identity 

Consider a meromorphic form TZ{x), with poles ri. 



f{x) = O (x"' ^) in all directions, and such that it has no singularities inside every 
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Then we have for x outside of all A-cycles (i.e. on the same side as ooq): 

n{x) = - ^ Res G{x, z)n{z)dz + ^ Vaix) (f n{z) dz. (4.61) 

proof: 

Since G{x, z) = l/{z — x) + . . . , we write Cauchy formula: 

n{x) = Res ^(x, z) n{z) dz (4.62) 

z— >x 

and we deform the contour of integration from a small circle around x, to contours 
enclosing all other singularities, i.e. the rj's and the Sj's. By doing so, G{x,z) has to 
cross the ^-cycles, and picks a discontinuity equal to 2i'jTVa{x) i.e. independent of z, 
so the contour integral of the product factorizes for each Aa- We thus arrive to: 

7^(x) = - V Res G{x, z)n{z)dz -T^TT- f z)n{z)dz 

i i •^'-^ 

+ ^Vaix)(b n{z)dz. (4.63) 
Then, we need to compute 

® G{x,z)lZ{z)dz. 
Jci 

Write that G{x, z) = ip^{z) K{x, z)/ip^{z), and integrate by parts: 

® G{x, z)TZ{z)dz = — f K{x, z)r{z)dz 

and using lemma 14.11 we see that this vanishes. 
□ 

5 Definition of correlators and free energies 

In this section, we define the quantum deformations of the correlation functions intro- 
duced in [T2l [T6] . Although the following definitions are inspired from (non hermitian) 
matrix models (see section [H]), they are valid in the present framework of an arbitrary 
Schrodinger equation, not necessarily linked to a matrix model. The special case of 
their application to matrix models will be discussed in section [91 

5.1 Definition of correlators 

Definition 5.1 We define the following functions Wn\xi, . . . ,Xn) called n-point 
correlation function of "genus" g by the recursio'^: 

Wf\x)=uj{x) , W!f\x^,X2) = B{X^,X2) (5.1) 



^here g is any given integer, it has nothing to do with the genus g of the spectral curve. 
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2g+2 

(9) 



(xo, J) = ^ X] r ^(^o,a;) (^iy„+2 (3^,3;, J) 

1 = 1 -^^i 



h=0 ICJ 

(5.2) 

where J is a collective notation for the variables J = {xi, . . . ,Xn}, and where ^ 
means that we exclude the terms [h = 0, 1 = ^) and [h = g,I = J), and where: 

W^:\x„ x.) = Wi^\x„ x„) - ^ (^7^ (5-3) 

i/ere xq and all the x'^s are outside of the A-cycles, i.e. on the same side as ooq- The 
contour Ci (defined in section 4-4-3\ ) is a contour which surrounds the branchpoint Li, 
i.e. a half-line of accumulation of zeroes, and chosen such that every sj is surrounded 
by exactly one Ci, and such that Ci doesn't intersect any A-cycle. Very often we shall 
write 

2g+2 

C=J2C^. (5.4) 

Appart from the precise definition of the kernel K, this definition is exactly the 
same topological recursion as in [T6j, a sum of residues around all branchpoints of the 
same expression. In other words, the topological recursion is independent of h. 

To shorten equation we will introduce the notation: 

To get: 
Theorem 5.1 

w!flixo, J) = ^j^d^ Kixo,x) Ui^\x, J) (5.6) 

proof: 

The only difference with the definition, is when we face a term like 

B{x,Xj)w!f\x, J/{j})- (note that there are twice this term). It can be split into two 

1 

(x — Xj)^ 

m TV., 

(x — Xj ) 



terms: B{x,Xj)Wn\x, J/{j}) and (^^J'^^-^2 Wn\x, J/{j})- The second term compen- 
sate exactly the d^^ ^"(}.^^^.)''^^ ■ Thus, the only difference between the two definitions 

-(9) 

(x — Xj) 

the same if these terms are null. This is the case because of Lemma 14. 1[ □ 



is the term: dx K{xo, x) dx^ u^l^ i^''^^ ■ Therefore the definitions are only 



35 



5.2 Properties of correlators 

The main reason of definition. 15.11 is because the wife's have many beautiful proper- 
ties, which generahze those of [12], and in particular they provide a solution of loop 
equations. We shall prove the following properties: 

Theorem 5.2 Each Wn^\xi, . . . with 2 — 2g — n < 0, is an analytical functions 
of all its arguments, with poles only when Xi Sj. Moreover, it vanishes at least as 
O (l/xf) when Xi ^ oo in all sectors. It has no discontinuity across A-cycles. 

proof: 

in appendix [B] □ 

Theorem 5.3 For all {n,g) ^ (0,0) we have 



We clearly have these properties for W2 (xo,Xi). By an easy recursion, the first 
property holds for xi, . . . , x„. The case of the variable Xq is special and requires expla- 
nation. Indeed for fixed values of xi, ...,Xn, the dependance in xq comes from K{xo, x). 
The theorem then comes from a permutation of integrals. Indeed, since the contour C 
never crosses any ^-cycles by prescription then we can permute the integrals in x and 
Xq. The nullity of the integral for K{xo,x) in 14.31 then gives the result. □ 

Theorem 5.4 For 2 — 2g — n < 0, the w!f^ 's satisfy the loop equation, i.e. Virasoro- 
like constraints. This means that the quantity: 




(5.7) 




(5.8) 



proof: 



2htiE^W^^l,{x, Xu Xn) + hd^W 




W':'{x,J/{3})-W':'{x,J/{j}) 



(5.9) 



is a polynomial in the variable x, of degree at most d — 2. 



proof: 



in appendix [Cl □ 
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Theorem 5.5 Each Wk^' is a symmetric function of all its arguments. 
proof: 

in appendix [El with the special case of W^^^ in appendix [Dl □ 
Theorem 5.6 The 3 point function W^"* can also he written: 

^(xi,X2,X3) = ^2^f^ yj^^ dx (5.10) 

(this can be seen as a quantum version of Rauch variational formula) 

proof: 

in appendix [D] □ 

Theorem 5.7 For 2 — 2g — n < 0, w!f\xi, . . . is homogeneous of degree2 — 2g — n: 




(5.11) 

proof: 

Under a change t^ — > Xtk, h — > \h, Cj — > Ae^, the Schrodinger equation remains 
unchanged, and thus ip is unchanged. The kernel K is changed to K/\ and nothing 
else is changed. By recursion, Wn"* is changed by A^~^^~". □ 

6 Deformations 

In this section, we will consider the variations of correlators Wn^ under infinitesimal 
variations of the Schrodinger potential U{x) or h. Infinitesimal variations of the re- 
solvent uj{x) can be decomposed on the basis of "meromorphic forms", and forms can 
be put in duality with cycles. The duality kernel pairing is the Bergman kernel. We 
will find in this section, that the classical ^ = formulae remain valid for /i 7^ 0, and 
generalize the corresponding form-cycle duality in special geometry. 

6.1 Variation of the resolvent 

Let's consider an infinitesimal polynomial variation: 
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where 5U is a polynomial of degree: deg5f/ < 2d. Since we have written U = V''^/4 — 
hV"/2 - P, we have: 

SU = —6V' - -5V" - —5V" - 6P (6.1) 
2 2 2 ^ ^ 

with 

d+l 



6V'{x) = J2^tkx''-\ (6.2) 

k=l 

and 6P is of degree at most d — 1: 

6P = {td+i5tQ + to6td+i)x'^~^ + lower degree. (6.3) 

Let us compute 6ip, or more precisely f = 6\nip = Sip/ip, let us write it: 

S^Ij{x) = f{x)tp{x) (6.4) 
The Schrodinger equation h^tp" = Uip implies: 

h\f4j)" - UftP = 6Uij - 5n^ i;" (6.5) 

i.e. 

n\f"i, + 2fi,') = {5U-2^-^U)^ (6.6) 



Multiplying by ip we get: 



h\fi,')' = {5U-2^-^U)i^' (6.7) 
n 



I.e.: 

<5(V7V') = f{x) = ^^-1^ i;\x') {SU{x') - 2 ^ U{x')) dx'. (6.8) 
therefore, since ui = V /2 + hip' /ip: 



2 ip{x) h%p\x) iooo ^ 



(6.9) 

If we write: ^ 

5^7 = — 5^' - -6V" - —V" - 6P (6.10) 
2 2 2 ^ ^ 

where (5P is of degree at most d — 1, and V /2 = uj — htp' /tp, we have by integration by 
parts: 



1 



X 



= rTHTZT / i^\x')(u{x')6V'{x')-6P{x') 

'ooo 



hip'^{x) 

6h{uj'{x')-l-V"{x'))) dx'. (6.11) 
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6.2 Decomposition of variations 

U{x) is a polynomial of degree 2d, it has 2d+l independent coefficients. If we assume 
that we have a solution of genus g < d — 1, this means that U is non generic, and 
satisfies d — 1 — g constraints. In the space of all possible f/'s, we shall consider the 
submanifold corresponding to U of genus g, which is a submanifold of dimension 

dim = rf + 2 + g (6.12) 

and we shall consider variations of U within that submanifold. Variations transverse 
to the genus g submanifold, are variations of higher genus and should be computed 
within a higher genus submanifold. 

Instead of the d + 2 + g independent coefficients of the polynomial U, it is more 
convenient to choose a system of "flat" coordinates in our genus g submanifold, given 
by: 

to,ti, . . . ,td+i, ei,...,eg. (6.13) 

We have indeed d + 2 + g coordinates. 

Let us write the variations as: 

d+i 

SU = J2 Ut, Stk + J2 + UhSh. (6.14) 

k=0 1=1 

6.2.1 Variations relatively to the filling fractions 

For the filling fraction 6ea we have 6V' = and thus: 

SU{x) = -6P{x) (6.15) 
where degSP < d — 2, so we decompose it on the basis of /i^'s: 

6P{x) = J2ca'K'. (6.16) 

and therefore, from eq. (16. 9p : 

5uj{x) = —''^^Ca' Va'{x). (6-17) 
a' 

Since 2m e a' = §a we have: 

2m5a,(x' = f StJ = — (p Ca"Va" = —Ca' (6.18) 
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This implies: 



and 



(6.19) 



6^^uli{x) = 2m Va{x) = <h B{x,z)dz. 

'Be 

(6.20) 

We shall say that the flat coordinate is dual to the holomorphic form Va-, which 
is itself dual to the cycle 3^- 



1 

2z7r 



UJ 



B. 



(6.21) 



6.2.2 Variations relatively to to 

We have: 

5U{x) = -6P{x) = -td+i x^-^ + Q{x) 
where degQ <d — 2. Using eq. fl6.9p we get: 

1 



bujix) 



i-t,+,x''-' + Q{x'))tP\x')dx' 



(6.22) 



(6.23) 



and the polynomial Q is chosen such that f^Su! = Oso that when decomposing Q{x) 
on the basis Va{x) and performing integrals over ^-cycles one finds the coefficients of 
the decomposition as integrals. Therefore we have: 



6uj{x) 



td+i Kd-i{x) 



(6.24) 



-td+l (j(d-l{x) Va{x) j kd^l{x')dx'- ^ Vo,{x) j ll)'\x')0i; 



where 



KJx) 



tjj^{x) 



x'^ V^^(x') dx\ 



(6.25) 



ooo 



and Kf^(x) is the A;"^ term in the large z expansion of K(x, z) = — Y2T=q ^Ik+t'^ com- 
puted in theorem 14. 2[ From theorem 14.51 we have G{x,oOi) = r]itd+i Kd-i{x). This 
shows that 



6toid{x) = G{x, ooo) = 77 {G{x, ooo) — G{x, oo_)) = / B{x, z) dz 



(6.26) 



where ooq is in the physical sheet, and oo_ is any infinity chosen in the second sheet. 



40 



We shall say that the flat coordinate to is dual to the 3rd kind meromorphic form 
—2G{x, ooq), which is itself dual to the chain [ooq, oo-]: 



tn = Resuj 



5taUJ = —2G{x, ooq) = I B{x, z) dz 

'ooo 



(6.27) 



where Res means the coefficient of in the given sector. 



6.2.3 Variation relatively to tk,k = 1 .. .d 



For k = 1, . . . , c? we have: 



V'(t\ 



degQ <d-2 



(6.28) 



and Q is chosen such that £^ 6uj = 0. Using eq. (16. 9p we write: 

6u{x) = (5a)(x) — '^^Va{x) (b 6uj{x') dx' 



where 



6lj{x) 



" V'ix' 



Since V'{x') = ^jtj+ia^'-', we have: 



x'^-' ij^x') dx' 



(6.29) 



(6.30) 



26uj{x) = tj+i Kk+j-1 

j=0 



(6.31) 



Let us compare it with the large z behaviour of G{x, z) in the physical sheet. We have: 

G(x, z) = V'(z) K{x, z) + 0{z-'^-^) (6.32) 



which means that the large z expansion of (^(x, 2) = "}2i],G]^[x) z ^ is given for k 



d by: 



and therefore 



j=0 



6uj{x) = --Gk{x) 

If we write the large z expansion of -B(x, z) in the physical sheet, we have 
5(x, z) = Y, Bk{x) z-^-^ = _ i ^ A; Gfc(x, z)z 



.33) 



(6.34) 



(6.35) 
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and thus 



1 z'' 
6t,uj(x) = — BAx) = Res — B(x, z) dz 
k ooo fc 



(6.36) 



We shall say that the flat coordinate tk is dual to the 2nd kind meromorphic form 
^ -Bfc(x), which is itself dual to a residue of B. 

6.2.4 Variations relatively to t^+i 

When k = d + 1, we have a few additional terms of degree > d — 2: 



Ut,,. (x) = ^ - ^ x'-' - tox'-' - Q{x) 



degQ<d-2 (6.37) 



and Q is chosen such that 6uj = 0. Using eq. 06. 9p we write: 

Suj{x) = 5uj{x) — ^^fQ,(x) ® Suj{x') dx 

where 



5Cj{x) 



1 



V ix ) 1^ dh id-i \ „i.2/i\ J I 



V{.X) ^ooo 

In other words we have: 



■ X 



X 



tox'"-^ ij\x')dx' 



(6.38) 



(6.39) 



2Suj{x) = ^tj+i Kd+j - dhKd-i - 'itoKd-i 

j=0 

Let us compare it with the large z behaviour of G{x, z). We have: 

2i 

G{x, z) = {V\z) -) K{x, z) - hd,K{x, z) + 0{z-''-^) 



(6.40) 



(6.41) 



which means that the large z expansion of G{x, z) = J2k ^k{x) z ^ is given for k = d+1 
by: 



Gd+iix) = - ^ tj+iKd+j + hdKd-i + 2toKd-i 

j=0 



and therefore 



6uj{x) = --Gd+i{x) 
If we write the large z expansion of B{x, z), we have 



B{x,z) = Y,B, 



k[X)Z ^ ^ = ~^J2^kix,z)z ^ 1 



(6.42) 



(6.43) 



(6.44) 
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and thus 

1 z'^ 
St uj{x) = —— - Bd+i{x) = Res — B{x, z) dz 

I (6.45) 

We shall say that the fiat coordinate td+i is dual to the 2nd kind meromorphic form 
5^+1 (x), which is itself dual to a residue of B. 

6.3 Variation relatively to h 

We have: 

5Mx) = - — ^ r ^\x') L'ix') - I V"{x') - 6nPix')) dx' (6.46) 
nip \Xj \ z / 

where 6nP is a polynomial of degree < d — 2 chosen so that £^ 6u! = 0. For the 
moment, we have not found a good way of writing this expression as an integral with 
B, and we leave that question for a future work. 

6.4 Form-cycle duality 

Notice that in all cases, except 5^, there exist a cycle Su* and a function A^^ such that: 



6u;{x)= / Bix,z)Al^iz)dz. (6.47) 

J Suj* 

We will use this generic notation later on in order to avoid specifying the 3 different 
cases. 

Under a suitable reparametrization z ^ z' such that dz' = A.}^{z) dz, we say that 
6u* in the variable z' is the cycle dual to the "meromorphic form" 6u!. 

6.5 Variation of higher correlators 

The following theorem allows to compute the infinitesimal variation of any Wn^^ under 
a variation of the Schrodinger equation. It tells about the "complex structure defor- 
mation" of our quantum Riemann surface. It can be regarded as special geometry 
relations. 

Theorem 6.1 Under an infinitesimal deformation U ^ U + 6U, we have: 

6Wjf\xu...,Xn)= I WiJi(xi,...,x„,x')A*(x')cia:' (6.48) 

J Suj* 

where {5uj* , A^^) is the dual cycle to the deformation of the resolvent uj ^ uj + 5uj. 
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proof: 

The loop equation for w!f\x, J) is: 

{2u{x) - V'{x) + hd,) W^'\x, J) + Ui^\x,x] J) = Pi'\x, J) (6.49) 

taking a variation 6 we have: 

{2(jj{x) - V'{x) + hd^) 5W^3\x, J) + {25u{x) - 5V\x))Wl^\x, J) + 5U!;^\x, x; J) 
= <5P(^)(x,J) (6.50) 

notice that 5Pn\x^ J) is a polynomial in x, of degree at most d — 2. 

On the other hand, consider the loop equation for wjf_l-^{x, J, x') and multiply it by 
A*(a;') and integrate x' along uj*, one gets: 

{2cu{x)-V'{x) + hd^) I w';,%{x,J,x')A*{x')dx' + I SUii^{x,x; J,x')A*{x')dx' 



(6.51) 



f Pl%{x,J,x')A*{x')dx' 

J UJ* 



That gives by recursion hypothesis for the computation of 
j^* 5t/^+i(x, x; J, x')h*{x')dx' and using [631 



{2uj{x) - V\x) + hd^) ( / l^iji(x, J, x')A*(x )rfx' - 5W^3\x, J) 



= 6P^s\x,J)- P!f^^{x,J,x')A*{x')dx' 

J U!* 

= Y,ai{J)K{x) (6.52) 

i 

where the right hand side is a polynomial of degee at most d — 2 in x, which can be 
decomposed on the basis hi{x). 

Solving the differential equation gives: 

/ I4^ii(x, J, x')K*{x')dx' - 5W^3\x, J) = Y^ ai{J) Vi{x) (6.53) 

% 

but since W^^ and lyiji are normalized on ^-cycles, this implies oti = 0, i.e.: 

iyi^\(x, J, x')A*{x')dx' = 6Wlf\x, J) (6.54) 

Jlj* 

□ 
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7 Free energies 

We use the variations and theorem 15.71 to define the Fg^s. 
Theorem 15.71 gives: 

/ d+l 

(^2-2g-n-hdn)Wi^^= to S*,, + 4 9,, + J] e.^,, ) W^f^ (7.1) 



fc=l i=l 



And in the previous section, we have seen how to write the derivatives of W^f '^ as 

n+l5 



integrals of WjfJi, that gives: 



^2-2g-n-hdn) Wjf^ = H. w!fl (7.2) 



n+1 

where H is the hnear operator acting as follows: 



H.f{x)=to / / + $^Res^/ + 5^6, i f. 



(7.3) 



Those equations allow to define H^q^^ = for n = and > 2 as: 
Definition 7.1 We define Fg for g > 2 such that: 

{2-2g-hdn)Fg = H.w['^ (7.4) 

It would remain to find the correct definitions of Fq (called the prepotential) and 
Fi. Fq and Fi should be such that under every deformation 5 = dt^, dt^, 9^. we should 
have 

6Fg = HsWi^\ (7.5) 

For example dF^/dtk = R.esx^uj{x)/k i.e. the coefficient of the term l/x^~^ in the 
expansion of uj{x) near ooq. 

We leave the definitions of Fq and Fi for a future work. 



8 Classical and quantum geometry: summary 

Let us summarize the comparison between classical algebraic geometry, and its quan- 
tum counterpart introduced here. 
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Figure 5: Classical case = of a two sheeted Riemann surface. The branchpoints 
are paired (in an arbitrary way) to form cuts, and the two sheets are glued along the 
cuts. Another possibility, is to draw a cut from each branchpoint to oo. The ^-cycles 
surround pairs of branchpoints in the physical sheet. There are also some degenerate 
branchpoints, which correspond to cuts of vanishing length. 



Summary 


classical h = 


quantum 




E{x,y) = J2i,jEijx'y^ 


^(a;, y) = Y.i,j Ei,jX' y^ , [y, x] = h 




E{x,y) = 


E{x, hd^)ip = 


n A/"nPTpl 1 1 "nf 1 r* a 1 
11 y ijci ciii jj ijn^cii 


= U{x) 


y'^ -U (x) , [y, x] = h, 


plane curve: 


deg U = 2d 




Potential: 


V\x) = 2(v/f/(x))+ 


2 sheets: 


y ±iV'{x) 


choice = ^0 \ in sector cxdq, ?7o = — 1 


resolvent: 


uj{x) = V'{x)/2 + y. 


u{x) = V\x)/2 + h^. 


physical sheet: 


y ~oo -^v'{x), U ~ to/x 


htp'/tp ~oo -iv'{x), U ~ to/x 

V 

sectors where ■?/; ~ e~2«' 


branchpoints: 


simple zeroes of U{x) 


half-lines of accumulations 




U{a,) = 0, U'iai) + 


of zeroes of 




z = l,...,20 + 2 


Li, 2 = 1,..., 20 + 2 


genus g: 


20 + 2 = 7^ branch points 






-1 < < d- 1 


double points: 


double zeroes of U [x) 


half-lines without accumulations 




u{h) = 0, u'idi) = 


of zeroes of 


genus g = — 1 


degenerate surface 


^^v/2h =polynomial 


^Q,-cycles 


surround pairs of 


surround pairs of half-lines 


a = 1, . . . ,g 


branchpoints 


of accumulating zeroes 


false ^Q,-cycles 


surround double 


links 2 sectors 


Q < a < d 


points 


where if) \ 


extra v4d-cycle 


surrounds last pair of 


surrounds last pair of half-lines 


a = d 


branchpoints 


of accumulating zeroes 


B-cjdes 




Ai n Bj = 6ij 


Holomorphic 
forms. 


= ^ 1 = TO C ^'(^') ^^(^') 

hi =polynomials, deg hi < d — 2 


1st kind 


normalized: §. Vi{x) dx = 6a,i, a = 1, . . . , 


differentials 


h^M = -i^a, ^/U'fk^) 


ip'^hi = 6a,i, a = 0+l,...,d-l 


Period matrix 







Summary 



classical h = 



quantum 



3rd kind form 



G{x, z) l/{z - x) 

z) = {2u{z) - V'{z) - hd,)K{x, z) 



Recursion kernel 



K(x.z) 



1 1 



K{x, z) 

K{X, Z) = K{X, Z) - Y,a Ca{z) 

C^{z) = §^^K{x',z)dx' 

V ' / htp'^(x) Jooo x'—z 



B{x, z) = — I dz G{x, z) 
B{x, z) ~ l/2{x- zf 



Bergman kernel 
2nd kind 



Symmetry: 



B[x, z) = B[z, x) 



§^^B{x,z)dx = 
B{x, z)dx = 2in Va{z) 



Meromorphic 
forms 



lZ{x)dx 



r\x)dx 



nx) = j^^Cr{x')^\x')dx' 



r{x) =rational with poles rj, r{x) = 0{x'^~'^] 
ReSr, r{x')i/j'^{x') = 



wiU^, J) = J:^i^ K{x, z) dz [Wi^-,'\z, z, J) 
where Cj surrounds the branchpoint Li 



Higher 
correlators 



W^^\x„(^l),X„{^2), 



Symmetry 



, . . . , 



a E Sn 



Variations and 
dual cycle 



U{x) U{x) + 5U{x) 
5U*: 5u{x) = j^,j,B{x,x') Ksu{x') dx' 



SV' = ^6tkX 



fc-i 



6t^uj{x) = ReSoo B{x, x') ^ dx' 



variation Stn 



K^i^) = Ij^ B{x,x')dx' 



variation (5e,- 



= B{x, x') dx' 



5W^'\x^, ...,xn) = J,^. • • - a:„, x') Asuix') dx' 



Variations of 
higher correlators 



T/T/ I ^ N _ X B{xi,z)B(x2,z)B(x3,z) 

yy3[Xl,X2,X3) — ^y,^^^ 



Ranch formula 



dz 



9 Application: Matrix models 

The reason why we introduced those Wn^^s is because they satisfy the loop equations 
for /3-random matrix ensembles. 

Consider a (possibly formal) matrix integral: 



Z=[ ciMe-^*^^^^) 



(9.1) 



where V{x) is some polynomial, and where E^^i = is the set of hermitian matrices 
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of size A^, is the set of real symmetric matrices of size and Ej\i^2 is the set of 

quaternion self dual matrices of size N (see |21j). 

Alternatively, we can integrate over the angular part and get an integral over eigen- 
values only f2T\ : 

Z= dX,...dXj, A(A)2/^ n e""^ ^^''^ (9.2) 

i=i 

where A(A) = ni<j('^i ^ ^i) Vandermonde determinant. 

This allows to generalize the matrix model to arbitrary values of (3. In particular, 
we shall choose /3 of the form: 



Notice that h = correspond to the hermitian case (3 = 1, and h —>■ —h corresponds 
to (3^ 1/(3. 

9.1 Correlators and loop equations 

Then we define the correlators: 



1 



Xk - A. 



(9.5) 



C 



and 

Wq = F = \^Z. (9.6) 

And we assume (this is automatically true if we are considering formal matrix 
integrals), that there is a large A^ expansion of the type (where we assume h, = 0(1)): 

oo 

iyfc(xi, . . . , xk) = J2iN/toy~''~''Wi'\xu ...,Xk) (9.7) 

5=0 

Wo = F = J2{N/tof-''Wl,'^ = J2Wtof-'^Fg. (9.8) 
a 9 
The loop equations are obtained by integration by parts, for example: 



48 



gives: 



Y 1 I o^>p 1 1 V'{\,) 

^\(x-Ai)2 j^^x-X,Xi-Xj to x-Xi 

y/ 1 I 1 1 V'iX.) 

■^^(x — Aj)2 X — Aj a; — Aj to x — Xi 

1-/3 , 1 1 iVv^ r(A,) 



^ / 1-/3 /?V — 



(a; — Aj)^ X — Aj a; — A,- tn x — Xi 

(,8 - ^)^W[(x) + ^(^WKx) + ^W^[x, x)) 



X - Xi 



(9.10) 

We define the polynomial 



AM - E ( "''1 (9.11) 



X — A 
We thus have: 

W^{x) + h— Wiix) + W2ix,x) = — {V'{x)Wi{x) - Pi(x)) (9.12) 

Using the expansion eq. fl9.7p . that gives the Ricatti equation 

iy^°^(x)2 + hd,wi^\x) = V'{x)W^''\x) - P^°^(x) (9.13) 
which is satisfied by uj{x): 

Wi'^\x) =uj{x). (9.14) 
generalizing to the integration by parts of 



k 

V^-^^M (9.15) 



I «i,...,ifej=i J j 

and using the expansion eq. fl9.7p to higher orders in to/N, one gets the loop equations 
of theorem 15. 4[ where 

Pk+i{x-x^,...,Xk) = ^(iV/to)'-'^-'Pifi(x;xi,...,Xfc) 



il,...,ik i=l ' c 
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In other words, the correlation functions of (3 matrix models, obey the topo- 
logical recursion of def. 15. 1[ 



Remark: 

In [15], a solution of loop equations for the /3-matrix ensemble was proposed, but 
that solution was such that U{x) was non-generic, corresponding to g = —1, and that 
had only a finite number of zeroes. This case implied that to was quantized. 
Generic matrix models cannot correspond to that situation. 

That solution was thus not very useful for actual matrix models. Here instead, we 
have the solution for every f/(x), i.e. every contour of integration for the Aj's, and 
therefore we have the solution of loop equations for the actual matrix model. 

9.2 Example: real eigenvalues 

Very often, we are interested in a matrix model with real potential V{x) of even degree 
(i.e. d is odd) and such that the eigenvalues are integrated along the real axis. The 
resolvent uj{x) is the Stieljes transform of the density of eigenvalues: 

Let us consider that it is defined by this integral in the upper half-plane for x G IHI+, 
and it is extended to the lower half-plane by analytical continuation. 

By definition, uj{x) is regular in the upper half-plane, therefore we look for a il){x) 
which has no zero in the upper half-plane, i.e. no zero on the half-lines Lq, Li, . . . , L^.i. 
I.e. it has at most d + 1 half-lines of zeroes , and thus: 



10 Non-oriented Ribbon graphs 

Consider the set of all closed connected ribbon graphs obtained by gluing the pieces 
represented in fig. [61 Closed means every half-edge is glued to another half-edge. 
Connected means every vertex is connected to any other vertex. See for example fig. [71 

We define the genus of such a ribbon graph Q as follows. We replace every twisted 
edge of ^ by a non-twisted one, we thus obtain another ribbon graph Q\ which is 
oriented. We define the genus of Q equal to that of Q': 

9{Q)=9{Q'). 
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The genus of Q' is computed as usual for oriented ribbon graphs, from the Euler 
characteristics of Q': 



xiG') = 2-2g = #{vertices(e?)} - #{edges(g?)} + #{single hnes(6;')} 



where single lines are the lines bordering each side of the ribbon edges. One should 
follow single lines and see how many connected single lines a graph contains. Obviously 
Q and Q' have the same number of fat vertices and fat edges (each edge containing two 
single lines), but they may have different number of single lines. 
This defines what we call the genus g of a ribbon graph. 

For a given Ribbon graph Q we call: 

• ni{Q) = ^unmarked vertices of degree i, for 3 < ? < (i + 1, 

• h{G) =size of the i^^ marked vertex, we have > 1. 

• e{g) = #edges, 

• q{G) = #twisted edges, 

• v{Q) = ^connected single lines, 

• 9{G) =genus, 

• ^Aut{Q) =symmetry factor of Q. 

Definition 10.1 Let 'n6^\v'), he the set of ribbon graphs Q with k marked vertices, q 
twisted edges, and of genus g, and such that Q' has v{Q') = v' connected single-lines, . 

Proposition 10.1 w6^\v') is a finite set. 



proof: 



The number of vertices of Q' is: 




The number of edges is twice the number of half-edges, i.e. 



k 




That gives: 



2 — 2g = ^{vertices} — #{edges} + v' = k 



i>3 i=l 



I.e. 




i>3 



1=1 
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unmarked 3 -vertex 



=j ^ unmarked 4 -vertex ^4 

^.'^ unmarked 5 -vertex ^5 
unmarked -vertex ty^ 



^■[^ marked vertex, size />0 x 

^ A K\ oriented, 

w A \ ^^^^ ^ marked edge 

= untwisted gluing 1/(3 

=^ twisted gluing 1-1/p 



Figure 6: Consider the set of ribbon graphs obtained by gluing those vertices. Marked 
vertices are of degree / > 1, they are oriented and have one marked half-edge. Un- 
marked vertices are unoriented, and are of degree > 3. Vertices are glued together by 
their half-edges, either twisted (with weight or untwisted (with weight 1 — 



Since the left hand side is fixed, we see that the number and size of vertices are bounded, 
so that there is only a finite number of possible oriented ribbon graphs Q'. Since Q' has 
a bounded number of edges, there is only a finite number of possibilities of twisting 
them, i.e. there are also only a finite number of graphs Q. □ 

10.1 Generating functions 

In order to enumerate the sets M.l^\v'), we define the following generating functions: 
Definition 10.2 We define: 

W^^\xi, . . . , Xfc; ts, . . . , td+i,(3; to) 

_ n-k/2 .v' ^ H ■■■ n-e{g) ( ^ _ ^\q{g) 

" k #Aut(^) . . . x^^^' ^ ^ 
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Torus Klein Bottle Projective plane 




Figure 7: Examples of ribbon graphs of genus (7 = 1. 



It is a formal series in powers of t^. 

Most often, for readability, we shall write only the dependence in the Xj 's; 

Wl'\xu . . . , Xfc; t3, • • • , td+uf3; to) = wl:'\xu • • • , a;^). 
Also, for k = we write 

10.2 Tutte's recursive equations 

Tutte's equation is a recursion on the number of edges to construct the ribbon graphs. It 
consists in finding a bijection between ribbon graphs of various ensembles, by recursion 
on the number of edges. Let M.\^^ be the set of ribbon graphs of genus g, and with 
k marked vertices of size li, . . . ,1^. 

Consider a ribbon graph Q G m[^^j^ ^ where L = {/i, . . . , Z^}, with marked vertices 
of degrees Iq + 1,L. 

Consider the marked edge of marked face 0. It is either twisted or untwisted. 
Several mutually exclusive situations may occur (see figE]): 

• on the other side of the marked edge, there is an unmarked vertex of size j + 1 
with j > 2. We then shrink the marked edge to concatenate the two vertices into one 
marked vertex of degree /q + j- The orientation is inherited from the initial marked 
vertex, and the marked edge is chosen as the first edge to the left of the shrinked edge. 
It is clear that we don't change the number of single lines in Q or Q'. We decrease the 
number of vertices and edges by 1, so we don't change the genus. We thus get a ribbon 
graph in M^^]^- j^, and this is weighted with weight t^+i (1 — = tj+i. 
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• on the other side of the marked edge, there is the marked vertex ? 7^ 0, of size 
k > 1. We then shrink the marked edge to concatenate the two vertices into one 
marked vertex of degree lo + li — 1. The orientation is inherited from the initial marked 
vertex, and the marked edge is chosen as the first edge to the left of the shrinked edge. 
It is also clear that we don't change the genus. Since we forget the marking of the 
other face, we shall get a symmetry factor /j, corresponding to the places where we 
glue to the i^^ marked vertex. We thus get a ribbon graph in M^^^j^ j^^^^ and this 
is weighted with weight k. 

• on the other side of the marked edge, there is the same marked vertex 0. Again 
we shall shrink the marked edge, i.e. shrink the 2 single lines. Several sub-situations 
may occur: 

* if the edge is untwisted, shrinking the 2 single lines splits the marked vertex of 
size lo + l into two vertices of size I' and Iq — I' — l- They inherit their orientation and 
marked edge from the initial marked vertex. We have increased the number of marked 
vertices by 1. The two new vertices are either connected together, or not. 

** If they are not connected, this means that the number of other marked vertices 
and the genus simply add up. We thus get two ribbon graphs in m/,^^' ^ ^i^Ziili l/l'i 
and this is weighted with weight 1/ (3. 

** If they are connected, we see that we get a new ribbon graph, with one more 
vertex, 1 less edge, and we have not changed the connectivity of single lines. The 
genus has thus decreased by 1. We thus get a ribbon graph in mI^^^\,_^ ^, and this is 
weighted with weight 1/(3. 

* if the edge is twisted, shrinking the 2 single lines doesn't split the marked vertex. 
Instead we get a new vertex of size 1^ — 1. We assign to it the orientation of the half- 
vertex situated left of the marked edge, and we mark the edge left of the initial one. 
We have decreased g by 1, and the genus is unchanged. We thus get a ribbon graph in 
M;5f2i,L) and this is weighted with weight (1 — 1//3) /q (indeed, there are /q places where 
we can glue the marked edge). 

For the generating function, those bijections read: 



d 





r{g-9') 

l+k-#L' 



,{x,Xl/l') 
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1 



VP 



pIU^.Xl) (10.2) 



we define 



d 

V\x) = ( X - ^tj+ix^' ) (10.3) 
i=2 



and the last term P^f^^^x] Xi) accounts for all tlie boundary terms, and it is necessarily 
equal to: 

Pif,(x; X,) = (V'{x) Wl:f,{x; X,))^. (10.4) 

This can be rewritten: 



V'{x)wtU^,X,) = Y^d^ 



1=1 



a',L' 

+wil-'\x,x,XL) 
+hd^wl:'^,{x,XL) 

+PiJi(a;,X^) (10.5) 

where 

^= V^-i/v^. 

In other words, the wjf^^s defined in section provide a solution to Tutte's equations. 
They are the generating functions counting our non-oriented ribbon graphs. One just 
needs to find the polynomial 

p(0) 

(x), i.e. U{x), and the choice of ip which is such that 
is a formal power series in to- 



ll Conclusion 

In this article, we have defined some "quantum" versions of quantities known in alge- 
braic geometry and applied them to the resolution of the loop equations in the arbitrary 
/9-random matrix model case, and in particular the enumeration of some non-orientable 
ribbon graphs. 

Our formalism recovers standard algebraic geometry and the invariants of [U] in 
the classical limit h 0. 

Instead of an albebraic equation, we have to deal with a differential equation, which 
we interpreted as a "quantum spectral curve", and we were able to generalize the 
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basic notions arising in classical algebraic geometry, like genus, sheets, branchpoints, 
meromorphic forms, of 1st kind, 2nd kind, 3rd kind, matrix of periods,... 

It is surprising to see that the notion of branchpoints become "blurred", a branch- 
point is no longer a point, but an asymptotic accumulation line. Also, there are two 
sheets, corresponding of the two possible large x asymptotic behaviors of tp{x) ~ 
exp ±V/2h, but in fact any solution is a linear combination of these two, so that we 
could say that we are always in a "linear superposition" of two states like in quantum 
mechanics. 

Another surprising thing, is that, in order for any cohomology theory to make 
sense, we need the cycle integrals of any forms to depend only on the homology class 
of the cycles, i.e. we need all forms to have vanishing residues at the Sj's. This "no- 
monodromy" condition is equivalent to a Bethe ansatz satisfied by the Sj's, like in the 
Gaudin model [2]. This provides a geometric interpretation of the Bethe ansatz, as the 
condition for cohomology to make sense. 

However, we still lack of a complete understanding of the situation, since most 
of our results explicitely depend on an initial sector 5*0 which we choose, whereas in 
algebraic geometry most of them only depend on the spectral curve and not on its 
parametrization. For instance the genus itself depends on a choice of sector. In some 
sense, the genus is no longer deterministic. 

Moreover, we still lack the proper definition of the spectral invariants Fg, indeed 
we have defined the F^'s only through solving a differential equation with respect to h, 
which is not as explicit as [IE] or [5]. Out of the -Fg's, we could expect the possibility 
to make the link with integrable systems and define a "quantum Tau-function" , like in 

m- 

Also, we restricted ourselves to the case of hyperelliptic curves, i.e. second order 
differential equations, or also a 1-matrix model. In a forthcoming paper, we shall 
generalize all this construction to arbitrary linear differential equations of any order, 
and generalize to a 2-matrix model. This work is underway, almost finished and the 
article is being written at this time. As for the hyperelliptical case, the notions of 
genus, sheets, branchpoints, forms, wjf^^s ... can be defined. Again there is a Bethe 
ansatz ensuring a no-monodromy condition so that all cycle integrals depend only on 
the homology class of cycles. So, there is no qualitative change, the difference is only 
technical, because the hyperelliptical case has big simplifications due to the involutive 
symmetry. The difference between the hyperelliptical case and the general case is 
comparable to the difference between [12] and [7], i.e. the definition of the kernel K is 
really more complicated, and there is a rather "big" technical step. 
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Then is would be interesting to see if the F^'s have some sort of symplectic in- 
variance, or more precisely some "canonical invariance", i.e. are unchanged under any 
change {x, y) (x, y) such that [y, x] = [y, x] = h. 

Finally, let us mention that we have developped a new geometrical approach to 
the study of D-modules, and it would be interesting to see how to relate it to more 
standard approaches, and also to the resurgence theory for studying the Schrodinger 
equation. 
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A Appendix: Proof of the loop equation for 

B{x, z) 

Let's first proove the first loop equation for z): Let's define: 

B{x, z) = \d. (2|^ - 9.) k{x, z) (A.l) 

i.e. we have: 

d-l 

B{x,z) = B{x,z) -^Va{x) (b B{x",z)dx" (A.2) 

1 JAr^ 



a=l 



Since (2^;^ + d^) Va{x) = ha{x) is a polynomial of degree < d — 2, it suffices to 
prove eq. fl4.35p for B{x,z). 
Let us compute: 



tjj{x) 2 ^(z) X — z 



2 ip(z)(x — z) (x — z) 

1 , g i^'iz) 



2 ' 



(x — z)^ ip{z){x — z) 
(A.3) 
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and therefore: 



(2^ + 4) (%,^)- ,, ,A +d.^ ^ = (A.4) 



This proves eq. fl4.35|] . with: 



pW(x,z) = / B{x",z)dx"- ^"(^) f dx" d^n^l)^{x")B{x",z). 

(A.5) 

Let's now proove the second loop equation for B{x,z): Similarly, let us compute 



{2§^ + d,)B{x,z 



Notice that the operator (j{z) = |(2^|;^ + d^) (2^^;^ — ^2), is equal to: 

U{z) = -\ dl + 2V{z)d, + \J\z) (A.7) 

which is also known in the litterature as the Gelfand-Dikii operator [8] (The Gelfand- 
Dikii differential polynomials Rk{U) are computed recursively by i?o = 1 and dzRk+i = 
U .Rk), which plays a key role in the KdV hierarchy. 

However, independently of any relationship with KdV, we get: 

IP - 1 

T^TT / '^^{x')dx' U{x').- 

V [x) J 000 ^ - z 

= ^ r ^^ix')dx' f- 3 + 

(A.8) 

We integrate the first term by parts three times, and we write Y = ip'/ip (we have 
Y' + = U): 

{2^ + dz)B{x,z) 

{x - zy ip-'ix) jooo - 



1 , Y(.) 1 r,^.(,,) 



(x — z)'^ {x — zY i/j'^{x) Jooo {x' — zY 
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J ooo M^;' - zf x' - z, 



t/j'{x')dx' — 



\2 



{x — zY {x — zY Ip'^ix) Jooo {x' ~ Z) 

tjj'^{x) V (x' — zY x' — z) 

1 , Y{x) Y\x) ^J_r ^2(^,) Y"{x') + 2Y{x')Y\x') 



(x — z)'^ {x — zY X ~ Z 1p'^{x) 7oog {x' — z) 

d Y{x) , 1 r , f2{U{z)-U{x')) , U'{z) + U{x') 



{x — zY dx X — z ip'^i^x) JooQ ^ (x' — zY x' — z 



(A.9) 
This implies that: 



'iIj{z) \ 2{x — zY ' dx x — z 

- -±- r ^^'(x')dx' I U'iz) + U{x'U 

~ r{x) Jooo ' \ {x'-zY x'-z ) 

= Q{z,x) (A.IO) 

which is clearly a polynomial in z. Taking integrals over x along Aa does not change 
its structure in z, and therefore: 

^'(-2) a N ^ 1 ^ ^ Y{x)-Y{z) 



(2 ^ I d^) (b{x,z) 

V'(^) ^ V ' 2{x — zY^ ' dx x — z 



ip'^{x) J^^ V {x' — zY x' — z 

= Pi'\z,x) (A.ll) 

is of the required form. 

By looking at the behavior of the various terms in the LHS of eq. fl4.36p when 
2; oo, we find that Pi'^\z,x) is a polynomial of degree at most d — 2 in z. 



B Appendix: Proof of theorem 15.2 



Theorem [O Each W^^\xi, . . . , x„) with 2 - 2g - n < 0, is an analytical functions 
of all its arguments, with poles only when Xi — > Sj . Moreover, it vanishes at least as 
O (l/xf) when Xi ^ oo in all sectors. It has no discontinuity across A-cycles. 

proof: 

We proceed by recursion on 2g + n. The theorem is true for W2^\ Assume it is 
true up to 2g + n, we shall prove it for Vr^^\(xo, Xi, . . . , 
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The integrand Un of theorem 15. II is singular only at x = s/s. As long as Xq is away 
from the s/s, we can continuously deform the ^-cycles and the contour C in order to 
have Xq outside of the ^-cycles, and the integral can be evaluated and is analytical in 
Xq. When Xq approaches Sj, we define Cj, a contour which surrounds all roots except 
Si, i.e: 

I = i +2iTT Res (B.l) 
Jc 

The integral over Cj can be evaluated and is convergent, thus it is analytical in xq. 

From the recursion hypothesis, all terms in the integrand are meromorphic in the 
vicinity of Sj, and thus the residue at Sj can be computed by taking a finite Taylor 
expansion of K{xq, x) = X]fc(^ ^ ^i)^ Ki^k{xo) in the vicinity of x ^ Sj. The result is a 
finite sum of terms of the type Ki^k{xo)- It is easy to see from the definition of K, that 
each Ki kixo) has only poles at Xq = Sj. Thus we have proved that has poles at 

the Sj's in its first variable. 

In the other variables, the result comes from an obvious recursion. 
□ 



C Appendix: Proof of theorem 15.4 

In this subsection we prove theorem 15. 4[ that all wjf^^s satisfy the loop equation. 

P!fl,{x,Xi...,Xn) = 2h^^W''^Ux,Xr...,Xn)+hdM'-^l,{x,Xu...,Xn) 

ICJ 



(CX) 



w'^:\x,j/{j})-wi'\x„j/{j} ) 

{X - Xj) 



is a polynomial in x of degree at most d — 2. 
proof: 

From the definition we have: 



J) = ^ (f dz K{x,z)Ui3\z,J) 
2nx Jc 

= ^<^^dzk{x,z)Ui^\z,J) 
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Then, notice that K{x,z) has a logarithmic cut along jooo;^:], and the discontinuity 
across that cut is: 

Ui^^ has no singularity outside of C, and thus we can deform the contour into a contour 
enclosing only the logarithmic cut of K{x^z), and therefore: 

^j^dz k{x, z) Ui^Hz, J) = -\ jy^^^ Ul?\z, J) 

We then apply the operator: + d^, that gives: 



D Appendix: Proof of theorem 15.6 



2m Jc Y'ix 

where Y[x) = — 2/i^^ 
proof: 



The definition of is: 



.(0) 

Ui vv 

AO), 



Wg ^(xo,Xi,X2) 

— — (p dx K{xo, x)B{x, xi)B{x, X2) 
«7r Jc 

— cp dx Kn G\ G'n 

4^7^ Jr 



(C.4) 



i2h^ + hd.,) l^j^dz kix, z) U^^^ (z, J) = -Ui^^ (x, J) (C.5) 
and therefore: 

Pii(x,J) = i2h^ + m,)w!flix,j) + ui^\x,j) 

= -^'^^ + ^^^^T."-^f/zC^iz)U^^\z,J) 

= - V h^{x) / dz C^{z) U!^\z, J) (C.6) 

which is indeed a polynomial of x of degree at most d — 2. 
□ 



Theorem 15.61 The 3 point function W^"* is symmetric and we have: 

„.{o)i ^ 4 /■ B{x,xi)B{x,X2)B{x,xz) 

'[xi,X2,x^) = — d) dx — — (D.l) 
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^ (( dxKo {{hK'l + YK[ + Y'K^){hK'^ + YK'^ + Y'K^)) 
Ait: Jc 

— (( dx Ko{ h^K'lK'' + hY{K[K'' + K';K'^ + hY'iK'K^ + K''Ki) 
+y2/s:;i^^ + YY\K^K'^ + K[K2) + Y'^K^K^ ) 



(D.2) 



where we have written for short Kp = K{xp,x), Gp = G{xp,x), and derivative are 
w.r.t. X. Note also that introducing Ki and K2 makes appear some additional and 
arbitrary logarithmic cuts from Xi to ooq and from X2 to ooq. But these cuts can be 
chosen arbitrarily since from the definition of W^^\xo, Xi, X2) it should not depend on 
that. Remember also that to use this definition of W^^^ we need to assume that xi 
and X2 are not circled by the contour C. Therefore we can choose the logarithmic cut 
of Ki and K2 inside the contour C like we have done it for xq. We now see that for 
example K0K1K2 has no singularity outside C and thus will not contribute because of 
theorem 14.11 Many other manipulations involving globally defined functions with no 
singularities outside C can be done. 

For example, using the Ricatti equation Y-^ = 2hY- + 4f/, we may replace Y-^ by 
2/iF/ and FiF/ by HYf. 

W^°\xo,Xi,X2) 

= — I dx Kq {hY{K[K'^ + K'lK'^) + hY\K'lK2 + K'^Ki) 
4^7^ Jc 

+2hY'K[K'^ + hY"{KiK'^ + K[K2) + Y'^K^K2 ) 

= i dx Ko ihY{K[K'J + hY'{KiK2)" + hY"{KiK2)' + F'^iTii^a ) 

4«vr Jc 

= -^j^dx Y'^KoK^K2 + h{Y"Ko{K,K2)' - {YK,)'K[K'^ - {Y' K,)'{KJ<2)') 
= dx Y"KoK,K2 - h{{YK,)'K[K'^ + Y' K'^{K^K2)') 

— I dx Y'^KoKiK2 - hYK'^K'Ji'^ - hY'{KoK[K'^ + KqKiK'^ + KKK2) 

AlTX Jc 

(D.3) 

This expression is clearly symmetric in xq^xi, X2 as claimed in theorem 15.51 

Let us give an alternative expression, in the form of the Verlinde or Krichever 
formula. 

^(0)/ ^ ^ N 2 / B{x,Xi)B{x,X2)B{x,x^) 



Wr{x^,x^,X2) = -(t dx ' ' ' ' (D.4) 

I'K Jc Y [x) 

proof: 

In order to prove formula ID.4[ compute: 

B{x,Xi) = -]f\x,x,) = -\g', = ^{hK': + YK[ + Y'Ki) (D.5) 
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thus: 



(D.6) 



1 
1 



dx 



B{x, Xi)B{x, X2)B{x, X3) 



dx 



Y'{x 



{hK'^ + YK'o + Y'Ko){hK'; + YK[ + Y'Ki){hK'^ + YK'^ + Y' K2) 



Y' 

+Y'^KoK^K2 



(D.7) 



1 / B(x,Xi)B(x,X2)B(x,X3) 

dx 



2m Tr Y'(x) 



-^f^dx hY{KoK[K'^ + KoK'{K'^ + K'.K^K'i + K'^K'lK2 + i^^i^i/^^ + KKK2) 
+hY\K'^K,K2 + iToir^'iTs + K^K.K'^ + ^iT^i^lir^ 

+^2(7^0^1^2 + KKiK'2 + KKK2) + YY'{K'oKiK2 + KoK[K2 + KqKiK'^) 
+Y'^KoKiK2 



Notice that = 2hY' + 4?7, thus we may replace Y^/Y' by 2^1^, and by 2^1^' and 
YY' by fiy", for the same reasons as before. Thus: 

— I dx ^(^'^i)^(^'^2)-B(a:,a;3) 



2i7r 7c r(a;) 

j dx hY{KoK[K'^ + KoK'lK'^ + K'^K^K'^ + K'^K'lK2 + KKiK'^ + K'^K[K2) 
+hY'{K'^KiK2 + i^o^i ^2 + KoK^K'i) + 2hYK'^K[K'^ 

+2hY'{KoK[K'^ + K'qK.K'^ + K'qK[K2) + hY"{K'QK^K2 + i^o^(^2 + K^K^K'^ 
+Y'^KqK^K2 

^jdx hY{Ko{K[K',y + K,{K'^K'^)' + K2{K'^K'J) 

+2hYK'^K[K'^ + Y'^KoKiK2 + h{Y'{K'^KiK2 + Koi^;fr2 + K^K^K'^))' 
Y^fdx hY{Ko{K[K'^y + Ki{K',K',y + K2{K',K[y) 
+2hYK'QK[K'^ + Y'^KoKiK2 

/ ^hYK'oK[K'2 + hY'{KoK[K'^ + i^^i^i^a + KKK2) 

lO'ZTr 
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Sin 

□ 



— i dx W^°\xo,xi,X2) (D.8) 
in Jc 



E Appendix: Proof of theorem 15.5 

Theorem E3] Each w!f^ is a symmetric function of all its arguments. 
proof: 

The special case of W^^^ is proved in appendix ID. II above. 

It is obvious from the definition that wjfJi{xo,Xi, . . . ,Xn) is symmetric in 
Xi, X2, . . . , Xn, and therefore we need to show that (for n > 1): 

wi'},{xo,x,, J) - l^ii(xi,Xo, J) = (E.l) 

where J = {x2, ■ ■ ■ , Xn}- We prove it by recursion on — x = 2g — 2 + n. 

Assume that every W^''^ with 2h + k-2<2g + n is symmetric. We have: 

iyii(xo,xi,j) 

= ^.Idx K{xo,x) (w!f-2'\x,x,xi,J) + 2 B{x, xi)W^'^\x, J) 

+2 EE <l,|(x,x„/)H^(!^,';)(x,j//)^ 



h=Q I&J 

(E.2) 



where means that we exclude the terms (/ = 0, = 0) and {I = J, h = g) . Notice 
also t 
have: 



also that wl^^2^ = wjf_^2^^ because n > 1. Then, using the recursion hypothesis, we 



w!fi,{xo,xuJ) 

2 (f dx K{xo,x) B{x,Xi)Wi3\x,J) 



+2 E E (^'' mtt'ii^i^'^ J /I) 
h I 

+2 E E ^i+i/i ^> mti\'^ J /I) 

h I 

+2 E E ^M/? ^/^) tel' ^) 



h leJ 
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+2 E E ///') 

h' I'd 

(E.3) 



Now, if we compute wjfJ^{xi,Xo, J), we get the same expression, with the order of 
integrations exchanged, i.e. we have to integrate x' before integrating x. Notice, by 
moving the integration contours, that: 

dx (p dx — (p dx (f dx = — (f dx — : Res (E.4) 



c Jc Jc 



2tii 



Moreover, the only terms which have a pole at a; = a;' are those containing B{x,x'). 
Therefore: 

Wif.ixo, xi, J) - xo, J) 

2 ^ dx iK{xo,x) B{x,xi) - K{xi,x) B{x,xo)) Wjf\x, J) 
—2 ® dx Res K{xo,x)K{xi,x') B{x,x') ( 

2W[Xn'\x\ X, J) + 2 E E ^Ih/? JIIW[%\{^\ I)) 

h leJ 

(E.5) 

The residue ReSjj/^a; can be computed: 

w!f_l,{xo,x^,J)-W!fl{x,,Xo,J) 
= 2 ^ dx {K{xo,x) B{x,xi) - K{xux) B{x,xo)) Wjf\x,J) 

- j) dx K{xo,x) — (^K{xi,x') 
2W[Xn'\x\x,J) + 2Y,Yl 

h l€J 

= 2 ^ dx iKixo,x) B{x,xi) - K{xi,x) B{x,xo)) J) 

— ^ dx K{xo, x)K'(xi, x] 



2wiii'\x,x,j)+2Y,Yl wtlJ^i^,J/m^1\Jli^J) 

h lej 

- j) dx K{xo,x)K{xi,x')—(^ 
2W[ll'\x',x,J) + 2Y,Y. 

2 j dx {K{xq,x) B{x,xi) - K{xi,x) B{x,xo)) Wi3\x,J) 
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— ^ dx K{xo,x)K'{xi,x) ^ 

2Will'\x,x,J)+2Y,Yl ^n-iiii^^J/nwl%^{x,I) 

h lej 

-- (j) dx K{xo,x)K{xi,x) — (^ 



h i&j 

(E.6) 

The last term can be integrated by parts, and we get: 
iyii(xo, xi, J) - Xo, J) 



2 (p dx {K{xo,x) B{x,xi) - K{xi,x) B{x,Xo)) Wjf\x,J) 

2Will'\x, X,J) + 2J2Y1 ^n-llli^^ J/mi^^ix, I) 



h leJ 

(E.7) 

Then we use theorem 15.41 

Wif^xo, xi, J) - Wii(xi, xo, J) 



= 2 ^ dx {K{xo,x) B{x,xi) - K{xi,x) B{x,xo)) W^^\x,J) 
+ j dx (^K'{xo,x)K{xux) - K{xo,x)K\xux)^ (^Pjf\x,J) 

+ {Y{x)-hd,)Wi'^\x, J) + ^a, 

(E.8) 



W!f\{x,,j/{x,}) 



X — Xj 



Since Pi^\x,J) and wjf}^{xj, J/{xj}) are entire functions of x, we can use the usual 
thorem 14.11 to say that they do not contribute. (Note again that we choose the loga- 
rithmic cut of Ki inside the contour C, and that we can do that because the contour C 
contains Xi.) 

wii(xo,xi,J)-iyii(xi,xo,J) 



2 (p dx (i^(xo,x) 5(x,xi) -fs:(xi,x) 5(x,Xo)) Wjf\x,J) 
{Y{x)-hd,)Wi'\x, J) 
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(E.9) 
Notice that: 



K'^K, - KoK[ = {GoKi - Kod) (E.IO) 



and B = —jG', therefore: 



w!f},{xo, xi, J) - Wii.ixu xo, J) 
-2 jf dx iKoG[ - K^G'o) Wjf\x, J) 

~2 jf dx (GoKi - J^oGi) iY{x) - hd,,)Wjf\x, J) 



(E.ll) 

we integrate the first hne by parts: 



(E.12) 

Notice that: 



wi%{xo,xuJ)-Wi'i,{xuXo,J) 
^ dx (K'qGi - K[Go) Wjf\x, J) 

+ j dx {KoGi - KiGo) W^3\x, jy 

jf dx (GoK, - i^oGi) - hd.,)Wj,3\x, J) 



K',G, - GoK[ = -J {KoG, - GoK,) (E.13) 

So we find 

w!f^,{xo,Xu J) - W^ii(xi,Xo, J) = (E.14) 
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Figure 8: When we shrink the single hnes of a marked edge, several possibilities may 
occur: 1) the other side is an umarked vertex of order j + 1, we get a new vertex of 
order Iq + j] 2) the other side is a marked vertex of order /j, we get a new vertex of 
order lo + li — 1; 3)4) the other side is the same vertex and the edge is untwisted. Then 
shrinking the edge splits the vertex into two vertices, this may disconnect the graph or 
not; 5) the other side is the same vertex and the edge is twisted. Then shrinking the 
edge doesn't disconnect the vertex. 
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